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I .   INTRODUCTION 

On  a  June  night  in  1948,  a  lone  ANSON  aircraft  -flew  a 
solitary  path  through  the  skies  o-f  East  Anglia,  England.   The 
wireless  operator  was  listening  -for  traces  o-f  signals 
propagating  through  the  dark  space  around  the  aircra-ft.   Near 
an  important  military  target,  a  signal  came  through  loud  and 
clear  in  his  earphones.   First,  a  series  of  dots  at  the  rate 
o-f  sixty  per  minute.   Then,  as  the  ANSON  droned  onward,  the 
dots  disappeared  into  a  steady  tone.   A  little  later,  the 
steady  tone  broke  up  into  a  series  o-f  dashes,  at  the  rate  o-f 
sixty  per  minute. 

This  signal  ,  known  as  a  Lorentz  beam,  was  transmitted  -from 
Germany  and  it  posed  a  considerable  threat  to  the  high  value 
military  target.   The  signal  was  really  composed  o-f  two  beams, 
one  beam  o-f  dots,  and  one  beam  of  dashes,  shown  in  Figure  1.1. 
By  using  the  reception  of  the  dots  and  dashes  as  an  indication 
of  whether  they  were  on  course,  German  bombers  could  navigate 
to  their  targets.  The  phasing  of  the  dots  and  dashes  was  such 
that  they  combined  into  one  steady  tone  in  the  overlap. 
Receipt  of  the  steady  tone  meant  the  aircraft  was  on  the 
course  laid  down  by  the  transmitting  site  and  receipt  of  dots 
or  dashes  meant  the  bomber  was  right  or  left  of  the  beam 
overlap.   With  that  information,  the  bomber  crews  could  easily 
get  back  on  course  and  perform  their  mission.  CI] 
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The  beams  were  ultimately  countered  by  the  British  using  a 
variety  of  methods,  but  the  invaluable  intelligence  gained  by 
the  ANSON  bomber  over  East  Anglia  was  only  the  beginning  of  a 
wartime  e-f-fort  that  was  to  save  thousands  o-f  lives,  both 
civilian  and  military.   That  brave  and  skillful  act  has  been 
repeated  hundreds  o-f  times  since,  in  a  "police  action"  and  an 
undeclared  war,  but  the  priceless  results  have  been  the  same, 
that  of  saving  lives. 


Kef.    wmz  Stesdj  note  zone 


nsmitter 


Dots  and  iishis  interlock 
so  that  where  beams  overlap 
a  steady  note  is  heard 


THE  LOHENZ    BEAM 


Figure  1.1.   The  Lorentz  Beam 

Since  the  arrival  of  the  first  electronic  transmitter  and 
receiver,  man  has  been  listening  to  his  own  radiations  for 
communications,  pleasure,  and  information.   World  War  II, 
however,  brought  a  new  maturity  to  listening  to  his  own 
radiations.   Now  he  was  listening  to  an  opponent's 
transmissions  to  derive  and  use  information  about  them. 
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About  1969,  the  term  Electronic  Support  Measures  (ESM)  was 
coined  C2:  p. 57]  and  now  reflects  the  important  business  cf 
obtaining  and  using  information  gained  on  the  electromagnetic 
radiations  of  another.  The  alert  operator  aboard  the  lone 
British  aircraft  used  his  early  and  undefined  ESM  receiver  to 
obtain  knowledge  that  was  used  against  the  German 
transmi  tters. 

A.   THE  ELECTRONIC  WARFARE  < EW)  ENVIRONMENT 

Of  all  the  sensors  that  an  enemy  might  use,  the  most 
important,  at  least  as  an  electronic  threat,  is  radar.  This 
sensor  is  almost  always  the  long  and  short  range  eyes  and  ears 
of  the  opposing  force,  and  denial  of  its  use  gives  a  strategic 
and  tactical  advantage  to  the  person  capable  of  doing  so.   E2: 
p. 7] 

Consider  the  case  of  a  penetrating  bomber.   In  order  to 
jam,  deceive,  or  otherwise  negate  the  effectiveness  of  the 
radars  of  the  defense  system,  emissions  must  be  detected, 
analyzed,  and  classified  as  to  their  threat.   In  this 
endeavor,  the  ESM  receivers  being  used  will  not  lack  for 
electromagnetic  sources  to  perform  their  analyses.   Early 
warning,  target  detection,  acquisition,  and  tracking,  or  even 
missile  and  interceptor  guidance  and  control  radars  will  paint 
the  aircraft  with  an  invisible  illumination. 

These  electromagnetic  energies  will  be  of  many 
frequencies,  antenna  scans,  pulse  shapes,  and  modulations,  and 
even  though  radar  systems  are  well  understood  and  sometimes 
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easily    countered,     the    sheer    numbers    and    types    o-f    received 
signals  will     tax    the    capabilities    o-f    the   ESM    system. 

As    the    bomber    nears    the    target,     the    signal     densities  will 
soar    because    the   defense    is    toughest    close    to    the    target.       For 
instance,    it    has   been    reported    that    over    188    SAM   sites  were 
once    located   on    the   Egyptian    side   o-f    the   Suez    Canal  .      As    each 
battery    transmits    several     signals,    easily   388-588    signals 
could   be   -found    -from    this    one   class   o-f    weapon.      When 
communication    and   radio    location    signals    are    added,    the    signal 
density   problem    is   even    -further    complicated.       C2:    pp. 27-33] 

Quite    easily,     the    combat    EW    environment    that    the    bomber    is 
in    will    reach    588,888    pulses    per    second.       C3:    p. 543 

Figure  1.2  C4]  gives  an  indication  o-f  the  pulse  densities 
the  bomber  might  encounter.  The  pulse  density  is  dependent 
upon  both  system  sensitivity  and  plat-form  altitude  but  is  not 
uniformly  dependent  on  these  parameters.  At  low  sensitivities 
<-28  to  -78  dBm)  pulse  density  varies  mostly  with  sensitivity 
and  not  with  altitude.  At  high  sensitivity  0-188  dBm)  pulse 
density    is    almost    totally    altitude   dependent. 

The    identification    o-f    various    emitters    and    the    evaluation 
o-f    the   possible    threats    they   pose    to    the    aircraft    requires    the 
sorting   of    these    pulses.      This   separating    or    deinter 1 eaving 
of    signals    is   a   difficult    problem    that    is   discussed    in    Chapter 
II  . 
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B.   ELECTRONIC  WARFARE  SUPPORT  MEASURES  <  ESM) 

The  broad  term  Electronic  Warfare  <EW),  or  sometimes 
re-ferred  to  as  Electromagnetic  War-fare  <  EMW)  ,  is  used  to 
denote  military  actions  that  involve  the  -friendly  exploitation 
o-f  the  electromagnetic  spectrum  and  the  prevention  of  hostile 
use  of  the  el ec tromaget i c  spectrum.   It  should  be  realized 
that  the  electromagnetic  spectrum  covers  from  zero  frequency 
to  the  infinite  frequency,  which  includes  the  optical, 
infrared,  and  laser  frequencies.   [5]   Coverage  of  these 
additions  to  the  RF  spectrum  adds  considerably  to  the  problem. 
EW  may  be  divided  into  three  major  divisions: 

1.  Electronic  Warfare  Support  Measures  (ESM) 

2.  Electronic  Counter  Measures  (ECM) 

3.  Electronic  Counter  Counter  Measures  (ECCM) 

Electronic  Support  Measures,  the  usual  shortened  form  of 
Electronic  Warfare  Support  Measures,  will  be  defined  shortly. 
ECM,  quite  briefly,  is  the  generation  of  intentional 
electronic  interference  between  electronic  systems  for  a 
useful  purpose.   ECCM  ,  also  quite  briefly,  is  action  taken  to 
counter  the  detrimental  effects  of  the  opponent's  ECM.  16: 
p.  13 

This  section  will  deal  only  with  the  ESM  arena. 

1 .   A  Definition  of  ESM 

Consider  a  specific  definition  of  ESM.   ESM  is: 

That  division  of  Electronics  Warfare  involving  actions  to 
search  for,  intercept,  locate,  record,  and  analyze  radiated 
electromagnetic  energy  for  the  purpose  of  exploiting  such 
radiations  in  support  of  military  operations.   C7:  p. 61 
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This  de-finition  carries  an  underlying  statement:  ESM 
is  the  source  o-f  information  required  to  carry  out  the  other 
Electronic  War-fare  divisions  as  well  as  threat  detection, 
warning,  avoidance,  target  acquisition,  and  homing.  C4: 
pp. 3-6]   In  addition,  by  examining  the  individual  actions 
specified  in  the  de-finition  we  can  gain  further  insight  into 
the  ESM  arena. 

The  search  tor  and  interception  o-f  electromagnetic 
radiations  is  simply  a  basic  part  o-f  the  ESM  mission.   The 
intercepted  signal  must  be  analyzed  and  a  determination  made 
o-f  its  character,  parameters,  and  location.   Locating  the 
emitter  is  usually  a  -function  -for  direction  -finding  antennas. 
Recording  o-f  the  signal,  certainly  the  parameters,  is  usually 
done  simultaneously  with  the  analysis.    On  some  missions  a 
recording  is  made  of  the  signal  for  postflight  analysis  with  a 
mimimum  of  information  being  displayed  to  the  operator  in  real 
time.   Generally  it  is  desired  to  produce  this  information 
-from  an  ESM  system  with  an  inherent  automatic  processing 
capability  that  detects,  classifies,  and  flashes  a  warning  if 
need  be.   C2:  p. 36] 

The  information  that  the  ESM  mission  provides  ranges 
from  that  collected  with  regularly  scheduled  and  dedicated 
collection  excursions  to  debriefs  of  crews  after  completion  of 
a  routine  patrol  mission.   To  a  tactical  commander,  this 
information  is  more  extensive  and  useful  than  plain 
intelligence.   It  provides  the  precise  state  of  the  electronic 
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de-fense,  including  technical  characteristics  and  emitter 
location,  and  becomes  part  of  mission  pre-f light  briefings. 

In  obvious  -fact,  radars  not  transmitting  are  radars 
that  can't  be  intercepted  and  analyzed.   Thus,  it"  there  are  no 
transmissions  -from  specific  desired  radars,  or  any  at  all, 
tactics  that  might  provoke  them  into  radiating  are  viable  in 
the  ESM  scenario.   Indeed,  this  action  is  included  in  the 
basic  rules  for  Electronic  Reconnaissance.   [2:  p. 65] 
2.   The  Role  of  ESM  in  Electronic  Warfare 

In  general ,  ESM  is  passive  electronic  warfare,  ECM  is 
active  electronic  warfare,  and  ECCM  may  be  either.   These 
distinctions  may  not  be  clear  cut,  however.   ESM  may  involve 
the  radiation  of  a  signal  to  determine  the  characteristics  or 
electronic  reaction  of  enemy  equipment.   ECM  may  involve  the 
passive  reception  of  enemy  signals  in  order  to  decide  which 
one  to  counter  C2:  p. 33.   The  latter  could  be  the  use  of  a 
tactical  ESM  system  integrated  with  an  ECM  system  for  platform 
protection.   The  relationships  between  ESM,  ECM,  and  ECCM  are 
shown  in  Figure  1.3.   C2:  p. 23 
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Figure    1.3.      The    Interactions    of    EW 
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The  previous  consideration  o-F  the  penetrating  bomber 
can  be  continued  in  the  context  of  ESM's  role  in  electronic 
war -fare. 

The  bomber  will  almost  surely  -face  an  integrated 
de-fense  network.   It  is  absolutely  essential  to  Know 
beforehand  the  technical  characteristics  and  location  o-f  the 
electronic  de-fense  that  will  be  used  against  the  bomber.   Some 
kind  o-f  electronic  reconnaissance  information  is  needed  in 
order  to  maximize  the  effectiveness  of  ECM  and  ECCM  equipment. 
Electronic  intelligence  (ELINT)  and  ESM  are  our  sources  of 
information . 

The  value  of  the  knowledge  gained  from  electronic 
reconnaissance  of  the  defense  system  is  proportional  to  its 
currency.   The  use  of  information  from  electronic  intelligence 
which  was  gained  from  many  sources  and  over  a  longer  period  of 
time  is  important,  but  the  immediate  needs  of  the  tactical 
commander  are  met  using  ESM.   This  information  is  collected  in 
the  hours  before  takeoff,  undergoes  minimal  analysis,  and  is 
used  for  mission  planning.   C2:  p. 37]   Thus,  ESM  is  electronic 
reconnaissance  to  determine  the  present  state  of  the  defense 
system. 

ELINT,  which  has  been  subjected  to  extensive  analysis, 
has  been  used  to  design  and  develop  the  ECM  suite  of  the 
bomber.   But  the  tactical  commander  still  uses  ESM  to  plan  and 
accomplish  the  mission.   Though  they  may  share  a  common 
collection  platform,  the  difference  between  them  can  be  seen 
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in  their  operational  control ,  end  use,  and  time  or  duration  of 
collection.   ESM  is  used  by  the  tactical  commander  today,  and 
ELINT  goes  to  the  rear  echelon  -for  further  study.   E2:  p. 74] 

Figure  1.4  [3]  shows  the  operational  role  of  ESM  in 
electronic  warfare.  Note  the  central  position  of  ESM  as  it 
"feeds"  both  ECM  and  ECCM  and  operational  readiness. 
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Figure  1.4.   Operational  Relationships  in  EW 
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C .   SUMMARY 

The  previous  sections  only  touch  upon  the  vast  nature  of 
EW  and  ESM.   They  provide  a  general  groundwork  for 
understanding  the  purpose  of  using  ESM  and  its  role  in  the 
electromagnetic  conflict. 

This  conflict  is  extremely  important.  New  and  modern 
weapon  systems  depend  upon  victory  in  the  electromagnetic 
conflict  as  a  prerequisite  for  victory  in  battle.   C2:  p. 13 

Underscoring  this  premise  are  new  policies  and  procedures 
by  the  Defense  Department  which  will  accelerate  development 
and  introduction  of  electronic  warfare  systems.   The  rapid 
escalation  and  introduction  of  new  Soviet  radars  places 
current  EW  sytems  at  a  disadvantage,  especially  in  countering 
fire  control  and  missile  guidance  radars.   [93 

The  future  probably  holds  a  greater  integration  of  ESM 
<certainly  tactical  ESM)  and  ECM  systems,  with  design  for 
military  aircraft  leading  the  way.   Already  systems  that 
combine  ECM  capability  with  passive  warning  equipment  for 
aircraft  are  being  sought  by  many  NATO  nations,  such  as  the 
British  Royal  Air  Force  C  18 3 ,  the  Belgian  Air  Force  [113,  as 
well  as  the  United  States  Air  Force  C 123 .   In  addition,  and 
with  considerable  impact  on  system  design,  a  trend  toward 
carrying  these  systems  internally  rather  than  in  pods  hung 
underneath  the  aircraft  seems  to  be  evident  C 18  and  123. 

In  the  hardware  area,  the  use  of  embedded  computers  and 
microprocessors  with  high  use  of  UHSIC  technology  will  become 
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prevalent,  along  with  all  of  the  ramifications  o-f  computer 
use.  The  sophisticated  signal  processing  needs  o-f  EW  will  have 
to  be  met,  with  verifiable,  reliable,  maintainable  so-ftware, 
and  the  development  o-f  high  capacity,  non-volatile  data 
storage .  C 13] 

Upcoming  and  -foreseeable  developments  are  EW  systems  -for 
space  applications  E2:  pp.  171-178].   All  of  the  unique 
environmental  conditions  o-f  space  will  generate  new  and 
di-f-ferent  challenges  -for  ECM  and  ECCM.   These  must  be  met  in 
order  -for  manned  space  stations  and  satellites  to  survive. 

D.   SCOPE  OF  THESIS 

A-fter  these  introductory  sections  Chapter  II  examines  some 
o-f  the  system  design  requirements  o-f  an  ESM  system.   The 
various  components  are  looked  at  in  some  detail .   Signal 
processing  is  addressed,  beginning  with  basic  parameter 
measurement,  and  ending  with  a  discussion  o-f  our  basic 
problem,  that  o-f  dein  ter  1  eaving  received  pulse  trains. 

The  attack  on  the  problem  begins  in  Chapter  III,  with  a 
discussion  o-f  the  theory  o-f  Walsh  and  Rademacher  Functions. 
Application  o-f  this  theory  is  carried  out  in  Chapter  Iv",  where 
the  Walsh  trans-form  is  used  on  the  dein  ter  1  eaving  problem. 
Finally,  conclusions  are  developed,  and  listings  o-f  the 
programs  used  are  given  in  order  that  follow  on  work  can  be 
done. 

The  attack  on  the  problem  using  the  Walsh  Transform  will 
begin  with  basic  examples  and  proceed  with  an  application  of 
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the  trans-form  directly  to  an  interleaved  pulse  train.  The 
goal  is  to  determine  whether  the  Walsh  Transform  outputs  a 
unique  feature  that  can  be  used  in  a  deinterl eavi ng  algorithm. 

The  radar  pulses  can  be  represented  by  time  of  arrival 
(TOA)  tags,  which  are  output  by  the  ESM  receiver  to  the 
preprocessor.  A  TOA  tag  can  be  thought  of  as  merely  a  "mark" 
in  time  that  indicates  the  receipt  of  a  pulse.   The  stream  of 
incoming  pulses  forms  a  string  of  TOA  TAG'S.   The  Transform  is 
applied  to  the  TOA  tag  string. 

The  Walsh  Transform  is  easily  applied  to  binary  or  two 
level  functions,  and  some  promise  is  held  for  its  use  in  a 
deinterl eaving  algorithm  if  a  PRI/PRF  recognition  feature 
exists  in  the  transform. 

With  the  speed  and  ease  of  computation  of  the  Walsh 
Transform,  a  processor  based  on  these  functions  should  be  fast 
and  simple.   Consider  the  programs  used  to  compute  the  Walsh 
transforms  in  this  thesis.   They  were  originally  written  in 
FORTRAN,  and  were  easily  converted  to  BASIC  for  use  in  this 
thesis. 

This  simple,  easy  to  use  language  normally  brings  a 
response  that  would  question  the  use  of  such  an  "elementary" 
language.   That  response  is  misinformed. 

BASIC  is  easy  to  learn,  and  easy  to  use.  It  performs  well 
when  compared  to  other  structured  languages,  and  has  excellent 
mathematical  and  scientific  processing  capability.  The  use  of 
the  language  certainly  has  not  hindered  the  growth  of  the 
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computing    power    available    in    a   microcomputer.       The    author's 
personal     computer,     the    IBM    Personal     Computer,    with    64K    o-f    RAM, 
has    the    same    amount    o-f    memory    that    a    -frontline    U.S.    Navy    P3-C 
Orion    ASW    aircra-ft    uses    to    process    seven    tactical    work 
stations,    plus    handling    all     the    navigation,    displays,    and 
other    systems. 

Throughout    this    thesis,    a   microcomputer    was    used    to 
compute    all     the   results    given,    and    only    the    lack    o-f    a    thesis 
quality    graphing    device    prevented    the    author    from    doing    all 
work    on    his    own    personal     computer.       Don't    underestimate    a 
microcomputer's    capabilities.      The    day    is    coming   when 
main-frame-like    computing    power    will     be    available    in    each 
o-f-fice    as    close    as    the    desk    or    table. 
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II .   ELECTRONIC  SUPPORT  MEASURES  SYSTEMS 

The  ESM  system  used  on  the  modern  battlefield  must  be  able 
to  give  rapid  and  accurate  assessments  o-f  the  complex  RF 
environment.    Detection  and  recognition  o-f  hostile  radar 
transmissions  must  take  only  a  -few  seconds  because  to  do  less 
invites  destruction  or  the  loss  o-f  the  ability  to  deal  a 
decisive  blow  to  the  enemy.   C3:  p. 54] 

The  design  o-f  ESM  systems,  in  view  o-f  their  importance, 
would  seem  to  -follow  a  rational  and  orderly  path  -from  concept 
to  development  to  operational  capability. 

This  usually  isn't  the  case,  un-for  tunatel  y ,  as  the 
reactive  nature  o-f  EW  wi  1  1  impede  an  orderly  -flow  or 
development.   Without  per-fect  knowledge  o-f  an  enemy's  radar 
and  EW  capabilities,  surprises  occur  which  tend  to  drive  the 
normal  development  into  a  "reaction"  driven  process.  The 
normal  process  also  tends  to  be  delayed  by  a  desire  to  wait 
until  a  precise  de-finition  o-f  the  threat  is  in  hand.   By  the 
time  it  is,  a  reaction  process  is  needed  to  counter  the 
threat.   C2:  p. 533 

As  re-ferred  to  in  the  Introduction,  page  25,  this  design 
process  is  undergoing  changes  to  improve  the  response  and 
lessen  the  reaction  design.   Changes  in  the  design  process 
will  not  change  the  design  requirements. 
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The  ESM  system  is  designed  to  intercept  many  different 
electromagnetic  signals.   Ideally,  the  system  should  be  able 
to  accomplish  as  quickly  as  possible  the  -following  tasks: 

1.  Intercept  a  transmitted  signal  at  any  frequency 

2.  Determine  the  type  of  modulation  in  the  signal 

3.  Identify  the  usable  intelligence  carried  in  the  signal 

4.  Measure  the  direction  of  arrival  of  the  waveform  for 
calculations  that  locate  the  transmitter. 

5.  Process  and  record  the  signal  characteristics 

6.  Display  to  the  operator  and  transfer  to  the  computer  for 
decision  making 

In  short,  the  system  must  gather,  process,  and  display  all 
signals  of  interest.   C2:  p. 1913 

Meeting  these  requirements  requires  a  lot  of  attention  and 
hardware.   No  single  tuner  can  cover  the  area  of  interest 
<usually  588  MHz  to  38  GHz),  so  many  tuners,  antennas,  and 
other  system  units  are  required. 

With  this  equipment,  an  ESM  system  approaches  the  problem 
with  these  basic  intents: 

1.  Intercept  the  signal.   Conversion  to  a  usable  form  must 
now  take  place  immediately.   An  alarm  should  be  given  to 
indicate  a  signal  has  been  intercepted.   Some  analysis 
or  sorting  must  be  done  immediately  to  determine  if  the 
signal  is  of  high  threat,  since  signal  frequency  and 
modulation  usually  correlate  with  degree  of  threat.   C2: 
p.  1911 

2.  Analyze  the  signal  for  further  information.   Once  it's 
decided  there  isn't  an  immediate  threat,  other 
parameters  of  the  signal  can  be  determined  for  recording 
and  later  decision  making. 

Each  of  these  approaches  has  more  or  less  importance 

depending  upon  the  platforms.   An  airborne  platform  has  less 

space  to  accommodate  equipment  for  reception  and  analysis,  but 

is  usually  exposed  to  more  immediate  threats.   Therefore,  an 
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airborne  plat-form  devotes  more  attention  to  interception  and 
immediate  analysis  than  the  ground  plat-form,  which  probably 
is  more  concerned  with  a  full  analysis  and  recording. 
Hope-fully,  this  entire  thesis  has  the  -flavor  o-f  applying  to  an 
airborne  plat-form. 

An  examination  o-f  a  typical  airborne  ESM  system  will  begin 
in  the  next  section. 

A.   COMPONENT  SYSTEMS 

Figure  2.1  C43  shows  a  typical  ESM  system.   O-f  particular 
interest  and  discussion  in  this  chapter  will  be  the  receiving 
system,  the  preprocessor,  and  the  processor.   The  Data  Files 
(Active  Emitter  and  Emitter  Parameters) ,  are  used  by  the 
system  as  records  -for  current  signals  being  received  (Active 
Emitter)  and  comparison  o-f  parameters  -from  processed  signals 
with  previously  entered  <a  priori)  parameters  o-f  signals  o-f 
interest  or  expected  signals  (Emitter  Parameters).   Electronic 
Counter-Measures  (ECM)  hando-f-f  and  equipment  is  not  really  a 
part  o-f  the  ESM  system,  but  is  included  as  a  -function  that 
could  be  available  -for  manual/automatic  jamming  or  deception 
o-f  high  threat  signals. 

1 .   The  Receiving  System 

Inseparable  in  their  functions,  the  antennas  and  the 
receiver  form  the  interception  and  detection  system  that  is  at 
the  forefront  of  the  entire  process.    This  section  will  cover 
a  few  of  the  important  aspects  of  the  antennas  and  receiver 
systems. 
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Figure    2.1.       An    ESM    System 
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a.   Antennas 

The  energy  received  by  any  plat-form  is  a  very 
small  -fraction  o-f  that  transmitted.   Consider  Table  2.1  C14: 
p. 9-33].   Notice  the  numbers  in  the  received  power  column  are 
magnitudes  below  one  milliwatt  in  almost  all  instances. 

Table  2.1.   Received  Power  Examples. 


Transmitter     Transmitter  Received 

Frequency     Peak  Power      Antenna  Gain     Range     Power 
(GHz)         (watts)  (dB)  (nm)      (dBm) 


5 

+  2 

288 

-38 

2888 

-58 

5 

-72 

18 

-24 

288 

-58 

i         187         40         iee      -7 

3  2  M  38 


18  1  38 

18  258  K  38 


48  188  K  35  18       -36 

High  sensitivity  in  the  system  is  needed  to  use 
this  power.   Good  antenna  design  is  the  starting  point  for 
reception  o-f  the  signals  that  are  desired  to  be  analyzed,  and 
good  design  is  \^ery   dependent  upon  the  application. 

The  basic  equation  that  is  -fundamental  to  the 
amount  o-f  power  available  for  use  is: 

Pp  ■  Pt6tx2/<<410  2R2)  (2.1) 

where  the  variables  and  units  are  as  indicated  below: 
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Pr  =  received  signal  power  in  watts 
Pj.  =  transmitted  power  in  watts 
G^  —  e-f-fective  transmitting  antenna  gain 
Gr  =  e-f-fective  receiving  antenna  gain 
X  =  wavelength  (same  units  as  R) 
R   =  separation  between  transmitter  and  receiver 

In  this  age  o-f  "if  you  can  see  it,  you  can  hit  it" 
accuracy  with  guided  missiles,  it  is  not  recommended  that 
shorter  ranges  be  used  to  improve  the  amount  o-f  power  being 
received.   Thus,  the  only  other  available  -factor  under  the 
designer's  control  (besides  receiver  design)  is  the  receiving 
antenna  gain. 

Our  antenna  gain  is  a-f-fected  by  several  -factors. 

G  =  4-SAP/X2  (2.2) 

where  A  and  p    are  de-fined 

A  =  capture  area 

P  =  antenna  e-f-ficiency 

The  capture  area  is  the  physical  size  o-f  the 
antenna.    The  e-f-ficiency  is  an  inherent  property  and  a 
■function  o-f  the  type  o-f  antenna  (usually  between  8.5  and 
8.6)  . 

The  beamwidth  o-f  the  pattern  that  has  the 
designer's  interest  also  has  an  e-f-fect  on  the  antenna  gain 
[15].   It  is  an  inverse  relationship   C 14:  p. 29-3]  that 
depends  upon  the  current  distribution  across  the  aperture. 

A  typical  re-flector  antenna,  -for  instance,  has 
G  %   28,888/(eB*B)  (2.3) 

where  9g  and  *g  are  the  ha  1 -f -power  beamwidths  in  degrees 
measured  in  the  principal  planes  o-f  the  pattern. 
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The    Gain    G    is    the    power    gain,    and   should   be    used 
in    radar    equations   because    it    includes    losses    introduced   by 
the   antenna    <Eq.    2.2).       Bearing    resolution    and    coverage 
■functions    -for    direction    -finding  would    lead    to    a    consideration 
of    the    directive    gain.      The   directive   gain    is  more    descriptive 
o-f    the   antenna   pattern. 

A    use-ful     expression    -for    directional     gain    is 
GD   ■    4*/<eB*B>  <2.4) 

where   Og    and   $g    are    the   half-power    beamwiths    in    radians. 
C  151 

Gain   must    be    considered    in    terms    of    signal 
reception.      Two   of    the    system's    tasks    are    to    receive    the 
signals  which   may    come   from    any    direction    (ornni 
considerations)    and    locate    the   emitter    (DF    considerations). 
The    two    tasks    are   not    necessarily   complementary,    so    a 
trade-off    in    system   performance    can    be    accepted    or    2    antennas 
can    be    used,    one    for    each    mission. 

For    an    airborne    platform,    weight,    space,    and 
structural    requirements  may    limit    the    size   of    an    ESM    antenna, 
but    the   altitude    of    the    aircraft   will    usually    offset    some   of 
these   disadvantages.       Dipoles,    slotted  monopoles,    and    surface 
wave   antennas    [ 14:    p. 29-38]    are    antennas    that    can    be   used    in 
airborne   applications. 

The   DF    function    of    the    system   could   use    any    of    the 
mentioned    types   of    antennas    if    the   system   uses    amplitude   or 
phase   comparisons   between    two    collectors.       C14:    p. 29-31] 
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Loop  antennas  are  a  common  amplitude  sensitive  DF 
antenna  o-ften  used  in  a  system  called  a  goniometer  C  16:  p. 
2h-143.   Horizontal  dipoles  or  dipoles  combined  with 
re-f lectors  are  also  antenna  types  that  are  used  in  DF  systems, 
as  well  as  Adcock  arrays  C 16:  p.2h-183 . 

Generally,  it  is  desired  to  have  high  gain  in 
antennas  -for  an  ESM  system,  and  the  speci-fic  application  will 
dictate  the  type  and  size  of  antenna  that  can  be  used.   More 
o-ften  than  not,  an  omni-direct ional  antenna  o-f  some  type  is 
used  -for  search,  and  a  rotating,  -fixed,  or  multi-port  highly 
directional  antenna  is  used  -for  direction  -finding, 
b.   Receivers 

Most  discussions  on  EW  wi 1 1  address  the  receiver 
technology  used  to  detect  the  signal  brought  in  by  the  antenna 
system.   Several  types  o-f  receivers  are  in  use  by  the  Air 
Force,  Navy,  and  Army,  with  the  mission  usually  dictating  the 
choice  of  technology. 

A  quick  and  general  overview  o-f  receiver  types  and 
technologies  -follows  in  this  section.   Although  some  o-f  the 
technologies  have  been  around  -for  some  time,  and  some  still 
yet  to  come,  they  are  increasingly  being  dominated  by  the 
processing  requirements  and  the  impact  o-f  microcomputers. 
Cost  saving  and  simplification  o-f  the  microwave  design 
requirements  are  a  result  o-f  this  increasing  influence.   In 
addition,  combinations  of  receiver  technologies  are  being  used 
to  overcome  limitations  of  individual  receiver  concepts. 
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The  -future  holds  new  technologies  -for  the  detection  of 
signals  and  increases  o-f  speed  and  bandwidth  o-f  single 
receiving  devices.   Surface  Acoustic  Wave  devices  will 
certainly  play  an  important  role,  as  well  as  acousto-optic 
techniques.   Increases  in  computational  speed  and  capacity  may 
aid  or  deter  some  technologies  by  eliminating  some  o-f  the  need 
-for  higher  sensitivity  in  the  detection  scheme. 

Table  2.2  C 17]  gives  some  typical  ESM  receiver 
spec  i-fications. 

Table  2.2.   Typical  Receiver  Specifications. 


Frequency  Range  8.5  to  18  GHz 

Signal  type  180  ns  pulse  to  CW 

Sensitivity  <-78  dBm 

Resolution  5  MHz  (pulse) 

1  MHz  (CW) 

Amplitude  accuracy  1  dB 

Bearing  accuracy  5*  RMS 

Pulse  width  resolution  188  ns 

Time  o-f  Arrival  resolution  58  ns 

Signal  density  18s  pps 

Intercept  probability  188X 


(1)  Direct  Detection.   Simplicity  in  design  gives 
the  direct  detection  receiver  advantages  in  cost,  reliability, 
size,  and  weight.  This  design  is  also  proven,  currently  in  use 
in  the  -front  end  o-f  radar  warning  receivers  (RWR)  ,  giving  a 
high  Probability  o-f  Intercept  (POD  -for  signals  above  its 
Minimum  Discernable  Signal. 

Today's  high  density  environments  pose 
problems  -for  the  DD  receivers,  degrading  their  performance 
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quickly.   Frequency  resolution  is  not  very  high  as  the 
bandwidth  is  usually  an  octave  or  more  in  width.   [18:  p. 26] 
Sensitivity  is  a  problem,  although  RF  amplification  can 
provide  improvements. 

Quite  susceptible  to  ECCM,  the  DD  receiver 
cannot  handle  PRI  agile  signals,  and  does  not  provide 
■frequency  measurement,  an  important  sorting  parameter. 
Figure  2.2  is  a  direct  detection  receiver. 


Threshold 
detector 


Processing 


Figure  2.2.   Direct  Detection  Receiver 


< 2)  Instantaneous  Frequency  Measurement  (IFM)  .   A 
■form  o-f  a  direct  detection  crystal  video  receiver,  the  IFM 
provides  -frequency  discrimination  and  measurement  with  a  high 
POI  .   O-ften  an  output  o-f  the  receiver  is  a  polar  display  o-f 
amplitude  versus  -frequency  C16:  p.2a-123  and  the  term  Polar 
Discriminator  is  used  in  place  o-f  IFM. 
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The  operation  of  the  discriminator  is  the 
heart  o-f  the  receiver.   Note  Figure  2.3(a).    Input  signals 
are  divided  into  two  parts  and  phase  correlated  after  one 
signal  is  delayed  by  a  known  time.  The  length  of  the  delay 
line  is  related  to  the  amount  o-f  -frequency  resolution  desired 
and  the  wavelength  that  corresponds  to  the  desired  bandwidth. 
By  banking  several  discriminators  and  different  lengths  of 
delay  line,  as  shown  in  Figure  2.3(b),  a  desired  frequency 
band  can  be  covered.   This  approach  allows  the  use  of 
practical  and  relatively  inexpensive  delays  rather  than 
attempting  to  use  one  delay  line  precise  enough  to  handle  the 
desired  frequency  region.   The  length  of  the  longest  delay 
becomes  the  frequency  resolution  of  the  bank.  C  18:  pp. 32-331 

Measurement  of  the  frequency  of  the  incoming 
signal  with  these  analog  discriminators  (or  current  technology 
Digital  IFM's)  allows  this  important  parameter  to  simplify  the 
sorting  process. 

On  the  down  side,  the  IFM  doesn't  handle  high 
density  environments  very  well  unless  a  subsystem  that 
determines  when  2  pulses  occur  simultaneously  is  included. 
Also,  the  delay  lines  require  special  attention  (constant 
temperature)  and  this  usually  increases  the  required  input 
power  to  the  system. 

Working  in  conjunction  with  an  analysis 
receiver,  the  IFM  makes  an  ideal  acquisition  receiver  by 
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Figure  2.3.   IFM  Receiver  Concepts 
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providing  the  -frequency  of  a  signal  o-f  interest.   Analysis 
receivers  o-f  another  type  can  then  be  tuned  to  the  frequency 
C18:  p. 32-34] 

(3)  Superheterodyne.   Most  common  of  the 
receivers  in  current  use,  the  superheterodyne  provides  high 
selectivity,  good  resistance  to  jamming,  and  is  a  proven 
design . 

Drawbacks  include  slow  time  in  scanning  the 
frequency  band  of  interest  and  inability  to  see  frequency 
agile  signals.   [19:  p. 4593 
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Figure  2.6.   Superheterodyne  Receiver. 

(4)  Channel ized.   Efforts  to  increase  the  POI  of 
superheterodyne  receivers  and  deal  with  an  inability  to  handle 
high  density  environments  led  to  the  concept  of  a  bank  of 
filters  followed  by  individual  signal  detectors  that 
determine  when  a  signal  is  within  the  filter  bandwidth.   The 
RF  spectrum  is  simply  broken  into  pieces.   See  Figure  2.5. 
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A  key  advantage  o-f  the  channelized  receiver  is 
removal  o-f  pulse  overlap  o-f  near  simultaneous  signals  <only 
pulse  frequencies  that  exceed  the  channel  spacing).   [17: 
p.  1853 

By  -following  an  initial  bank  o-f  -filters  with 
multiplexers  and  -fixed  oscillators,  then  repeating  the 
scheme  with  a  smaller  bandwidth  -filter,  the  RF  input  can  be 
down  converted  to  a  -final  baseband  -frequency  range  with  a 
higher  -frequency  resolution. 
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Figure  2.5.   Channelized  Receiver 

<5)  Compressive.   This  receiver  can  be 
characterized  as  a  -fast  scanning  superheterodyne  receiver. 
Ideally,  the  local  oscillator  scans  the  RF  bandwidth  being 
covered  in  a  time  less  than  the  narrowest  pulse  to  be 
intercepted.   Its  high  POI  and  the  ability  to  handle  wideband 
signals  and  -frequency  agile  signals  makes  this  receiver  a  good 
choice  -for  an  ESM  receiving  system.  It  has  an  excellent 
ability  to  separate  signals  closely  spaced  in  frequency. 
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Basically,    the   compressive    receiver    provides 
an    IF   signal     that    is   up-chirped    in    frequency    -from    the   RF 
input.      The   compressive   filter    operates    in    such    a  way    that    the 
low   frequency    components   of    the    IF    signal     are   delayed    longer 
than    the   high    frequency    components,    with    the    result    that    the 
shape   of    the    output    pulse   of    the    detector    is    in    essence    the 
Fourier   Transform   of    the    IF    signal     pulse.      The  maximum 
amplitude   of    the    detector    output    falls    at    a   point    in    the    scan 
time    that    is   proportional     to    the   frequency   of    the   RF    input. 
[18:    pp.    38-31] 

Critical i ty   of    the    alignment    between    sweep    and 
compressive    delay  makes    this    complex    technique    difficult    to 
manage    C 19:    p.    459].       But    it    is    at    the   forefront    of    some    of 
the   newest    technologies,    with    SAW's   performing    the   delay 
necessary    for    the   up    and    down    chirp.      Note    Figure   2.6. 
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Figure  2.6.   Compressive  Receiver. 

<6)  Bragg  Cel 1 .   An  acousto-opt ic  Bragg  cell  will 
interact  a  propagating  acoustic  wave  with  an  optical  beam  to 
give  a  diffracted  output  proportional  to  the  frequencies  and 
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power  present  in  the  RF  input  signal.   The  acoustic  wave, 
which  varies  the  index  o-f  re-fraction  of  the  cell  material  ,  is 
generated  by  applying  the  RF  input  to  an  acoustic  transducer 
in  the  Bragg  cell.   A  laser  beam  is  de-flected  and  modulated  in 
intensity  by  the  index  variations  which  are  proportional  to 
the  RF  input  -frequency.    A  Fourier  Trans-form  lens  collects 
the  proportionately  di-f-fracted  light  which  then  -falls  on  a 
photodetector  array.   The  output  o-f  the  detector  array  is 
similarly  proportional  to  the  RF  input  -frequency.   C  18:  p. 283 
Having  pluses  -for  high  sensitivity,  high  POI  ,  and 
high  selectivity  makes  this  receiver  a  choice  o-f  the  -future. 
It  doesn't  handle  -frequency  agile  sources  well,  however, 
because  o-f  the  time  needed  in  determining  the  output  o-f  the 
photodetector  array.   Pulse  width  measurement  is  also 
eliminated  with  the  use  o-f  the  detector  array  because  the 
array  only  responds  to  the  energy  contained  in  several  pulses 
■from  the  lens.   C19:  p. 4683   Figure  2.7  is  the  receiver. 


V 

Xfm- 
Lens 

Sensor 
Array 

* 

i 

Laser) 

— ► 

Channeli zer 

— ► 

Filter 

— ► 

Bragg 

te 


Processing 


Figure   2.7.      The   Bragg   Cell    Receiver 
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2.   The  Processing  System 

A  steady  stream  o-f  pulse  descriptor  words  -from  the 
receiving  system  must  be  analyzed  -for  the  vital  in-formation 
that  is  contained  in  them.   The  processing  system  must  extract 
this  in-formation  and  make  decisions  regarding  the  emitters 
whose  signals  it  has  processed.   In-formation  extraction, 
analysis  and  identification  must  be  -fast  to  avoid  a  mission 
and  li-fe  threatening  situation. 

This  section  will  discuss  the  overall  purpose  and 
general  operation  o-f  preprocessors  and  processors  in  the  ESM 
system.   A  more  complete  discussion  will  be  -found  in  Re-ference 
28. 

a.   Preprocessor 

Basically,  the  preprocessor  must  prepare  the 
signal  received  -from  the  receiving  system  -for  prel  iminary 
analysis  and  -for  advanced  analysis  in  the  main  processor. 
Figure  2.8  is  an  ESM  preprocessor.   C28:  p. 83 

A-fter  digitization,  the  initial  parameters  o-f 
pulse  width,  time  o-f  arrival,  direction  o-f  arrival,  amplitude, 
and  -frequency  must  be  measured  on  a  pulse  by  pulse  basis. 
(See  Parameter  Measurements  in  this  chapter.) 

The  pulse  data  must  be  compared  with  previously 
received  and  analyzed  signals  to  see  i-f  this  signal  has 
already  been  processed.   Pulse  data  is  compared  with  signal 
data  in  the  -files  that  contain  the  active  emitters  and  the 
uninteresting  emitters.    I-f  a  match  is  present,  then  the  data 


45 


V 

c 

K 

UJ 

Ul 

3 

t- 

G 

s 

Ui 

C 

a 

c 

IT 

a. 

c 

a 

r 

X 

►«■ 

UJ 

c 

01 

►» 

a. 

K 

<T 

Ul 

U 

M 

_i 

^ 

UJ 

c 

|M 

w 

3 

4 

a. 

13 

r 

•* 

a. 

0) 

2 

j 

u. 

r 

j 

L 

e 

u 

or 

> 

o 

c 

UJ 

Ul 

U. 

in 

J 

ac 

UJ 

K 

a. 

(j 

ui 

a 

*- 

a. 

2 

a. 

U 

a 

UJ  c 

>-  t-         Ul         3 

Ul  —          _l           _ 

E  £       —       X 

«  —       u.       u 

a:  -i 

i  a 

a.  — 


►-  —  •-  _i 

Z  *-  t-  — 

—  ji  —  u. 

Z  Ul  £ 


A 

i 

c 

a      r 

UJ 

O        uJ 

M 

-J        u. 

■. 

«      u. 

►- 

Z        3 

-* 

<t       a 

a 
a 

i 

i 

z 

IT 

o 

> 

X 

u 

U. 

il 

Figure    2.8.       Preprocessor 
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can    be   deleted    -from    the    process.       Reduction    o-f    the    data    stream 
and   sorting   of    the   signals    is    the   purpose   and   result    o-f    this 
deletion    in    addition    to    avoiding   unnecessary    processing.       One 
o-f    the  main    jobs   o-f    the   preprocessors    is    to   convert    the 
approximately    18s    to    186   pulses   per    second   being    received    to 
about    1888    pulses   per    second    that    a   good  main    processor    can 
handle.       [21:    p. 164] 

Unmatched    or    new    signals  must    be    deinterl eaved 
(see   deinterl eaving    in    this    chapter)    and   processed    in    the   PRI 
processor    section.       Calculation    o-f    the   PRI     is   usually    done 
simply   by   subtracting  TOA's    o-f    similar    pulses.       Staggered 
PRI's   and    jittered    PRI's   have    to   be   handled  with    a  more 
complex    algorithm.       [28:    p.     12] 

With    a   known    PRI,    TOA   prediction,    additional 
sorting,    and    comparison    o-f    calculated   PRI's    is    done   until     a 
PRI    is   con-firmed    to    exist    in    the    data.    This    in-f  ormat  ion ,    the 
newly   calculated   PRI    and    the    initial    parameters   o-f    the   pulse, 
is    then    passed    to  main   memory    o-f    the   processor    -for    analysis 
and   updating   o-f    the    active   -files   so    that    pulses    o-f    the 
incoming   stream   can    be    deleted   as    soon    as   possible. 

In    summary,    the   preprocessor    prepares    the    incoming 
signals   -for    processing,    measures    their    initial    parameters, 
per-forms   preliminary    analysis   and   sorting,    and   passes    the    data 
information    to    the  main    processor. 
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Figure  2.9.   Main  Processor 


48 


b.   Main  Processor   C29:  pp. 13-183 

A-f ter  al  1  the  action  in  the  preprocessor  and  the 
data  reduction  that  has  been  accomplished,  there  is  still 
plenty  to  do.   The  -functions  o-f  the  main  processor  are: 

1.  Advanced  analysis  o-f  the  signal  -for  scan  rate  and  exotic 
emitter  (such  as  chirp  radars)  identification 

2.  Emitter  and  plat-form  iden  ti-f  ication 

3.  Updating  the  emitter  -files 

4.  System  control  and  monitoring 

5.  Inter-facing  with  the  operator  and  ECM  equipment 

Since  the  processor  is  usually  a  general  purpose 
computer,  additional  -functions  o-f  navigation  or  routine  system 
calculations  could  be  added.   It  can  also  help  the  PRI 
processor  handle  its  calculations. 

Emitter  iden  ti-f  ication  is  accomplished  by 
comparing  the  measured  emitter  parameters  with  emitter 
parameter  lists  until  an  identification  is  made  or  the  signal 
is  tagged  as  unidentified.   Interface  with  the  operator  for 
manual  analysis  is  often  the  result  of  an  unidentified 
signal . 

3.   Display  System 

An  operator  must  have  a  visual  presentation  of  the 
signal  and  its  parameters.   Almost  universally  CRT  displays 
are  used.   CRT's  usually  have  an  intensity  modulated  or  a 
deflection  modulated  design,  with  advantages  and  disadvantages 
in  each  one. 

Multi-gun  CRT  displays  can  create  a  combined  display 
that  offers  a  pulse  analyzer,  presenting  frequency,  PRF,  PW, 
and  the  presentation  of  an   amplitude  spectrum. 
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A  warning  system  usually  operates  in  conjuction  with 
detection  o-f  high  threat  signals.   Flashing  lights  and/or  an 
aural  alarm  will  occur,  indicating  immediate  attention  is 
required. 

The  audio  system  operates  in  parallel  with  the  warning 
system.   Offering  not  only  alarms  to  high  threat  signal 
reception,  the  system  usually  has  accomodations  that  allow  the 
operator  to  monitor  the  received  signal  through  headsets  or  a 
speaker.   A  trained  and  experienced  operator  can  determine 
scan  rate  and  analyze  the  scan  modulation  type  i-f  the  signal 
density  is  low  enough. 

Bearings  from  DF  equipment  are  usually  presented  on  a 
polar  display.   The  direction  is  usually  indicated  in  relative 
bearings  o-f  8  to  368*  around  the  plat-form,  with  some  output  or 
means  o-f  determining  the  true  bearing  for  plotting  purposes. 

Computer  management  of  an  ESM  system  allows  for 
most  calculations  of  bearing,  PRF ,  PW,  and  frequency  to  be 
accomplished  and  shown  as  soon  as  the  bearing  is  displayed  to 
the  operator.   An  excellent  example  of  computer  managed 
displays  is  the  ALQ-78  ESM  system  on  the  U.S.  Navy  P3-C 
Orion.    Digital  readouts  of  the  PRF,  PW,  and  frequency  are 
available  next  to  the  displayed  bearing  from  the  aircraft. 
The  bearing  remains  on  the  screen  and  updates  automatically 
until  removed  by  the  operator  or  the  system.   With  several 
cross  bearings,  the  emitter  location  can  be  obtained  since  the 
computer  is  not  only  managing  the  display  but  the  navigation 
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of  the  aircraft.  By  "hooking"  the  place  where  the  bearings 
cross,  the  operator  can  cause  the  computer  to  calculate  the 
lat/long  of  the  emitter  fix. 

B.   ESM  SIGNAL  PROCESSING 

The  acquisition  of  signals  is  usually  a  design  problem 
that  is  driven  by  the  characteristics  of  the  signals  to  be 
acquired.   The  designer  chooses  antennas  and  receivers  to 
detect  signals  with  regard  to  the  mission,  the  platform,  and 
the  frequency  band  of  interest.   Generally,  he  desires  the  ESM 
system  should  have  rapid  acquisition,  high  sensitivity  and 
frequency  resolution,  and  wide  frequency  coverage. 

After  these  signals  are  acquired,  signal  processing  is  the 
order  of  the  day.   Processing  is  required  immediately  to 
determine  the  high-threat  signals  <especially  for  tactical 
aircraft) ,  and  a  lot  of  processing  is  required  in  today's 
dense  electromagnetic  environments.   The  signals  are  getting 
more  complex  also,  with  jittered  and  staggered  PRF's  and 
variable  PW's  in  addition  to  frequency  agility  and  new 
schemes  of  radiating  and  detecting  a  return. 

Computational  speed  becomes  one  answer  to  solve  these 
processing  riddles,  and  the  debate  between  digital  and  analog 
techniques  fuels  upon  this  quest  for  speed. 

Surface  Acoustic  Wave  <SAW)  devices  and  Charge  Coupled 
devices  (CCD's)  are  being  challenged  by  digital  technology 
staked  by  the  Pentagon's  Very  High  Speed  Integrated  Circuits 
Program  <VHSIC> .   These  digital  IC's  are  being  structured  to 
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per-form  high  speed  calculation,  with  a  N  X  N  matrix  inversion 
<a  N3  operation)  being  a  gauge  o-f  success-ful  application. 
Today's  analog  technologies  currently  o-ffer  high  speed  also, 
but  will  require  more  precision  and  dynamic  range  to  meet  this 
challenge.   [22:  p. 91] 

In  the  end,  the  ESM  system  designer  will  win  as  new 
digital  architectures  and  new  device  structures  will  yield 
increase  in  speed  and  throughput. 

This  section  will  discuss  important  signal  processing 
parameters  o-f  the  ESM  system,  and  the  techniques  used  to 
separate  the  signals  -for  identification  and  processing  -from 
the  immense  amounts  o-f  pulses  being  received.   The  entire 
process  is  a  continuous  -flow  o-f  calculations,  comparisons, 
and  decisions  involving  tolerances  in  -frequencies,  pulse 
width's  <PW)  ,  pulse  repetition  intervals  (PRI) ,  and  other 
parameters . 

1 .   Parameter  Measurement 

Listed  below  are  the  initial  parameters,  as  provided 
by  the  receiving  system  to  the  preprocessor.   They  are  the 
basic  tools  that  are  used  to  deinter leave  (see  next  section) 
the  pulses  and  match  them  to  existing  tracked  or  not  tracked 
emi  t ters. 

1.  Direction  o-f  Arrival 

2.  Ampl i  tude 

3.  Frequency 

4.  Pulse  Width 

5.  Time  o-f  Arrival 
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Those  that  aren't  matched  must  be  analyzed  -further  by 
calculating  a  pulse  PRI  and/or  undergoing  advanced  analysis. 
Calculation  o-f  the  PRI  is  done  in  the  preprocessor. 
Advanced  analysis  is  accomplished  in  the  main  processor  using 
derived  higher  order  parameters  during  the  task. 

A  discussion  o-f  some  aspects  o-f  how  the  initial 
parameters  and  higher  order  parameters  are  measured  -follows  in 
this  section.  [21:  p.  1663 

a.  Direction  of  Arrival  <DOA) 

An  important  parameter,  the  DOA  could  be  measured 
with  a  rotating  DF  antenna  but  this  greatly  lowers  the 
Probability  o-f  Intercept  <POI)  .   A  more  usual  method  uses  a 
multi-port  (antennas)  bearing  measurement  subsystem  that  has  a 
high  <*1887.)  POI  . 

Phase  comparison  techniques  can  also  be  used, 
providing  greater  accuracy  but  with  a  tradeo-f-f  o-f  smaller 
bandwidth . 

b.  Ampl i  tude 

Usually    the    amplitude    is    termed   by    the   peak   value 
o-f    the   received   pulse  with    some    tolerances    included    to   provide 
for   variances    in    emitter    radiated   power. 

The   amplitude    parameter    is    used    in    the    advanced 
analysis    section    of    the    typical    processor    to   determine    a   scan 
pattern,    but    is   not    considered    a    good    sorting   parameter 
because   of    propagation    corruption,    multi-path,    and   other 
variables    that    influence    the   amplitude   of    the   received   pulse. 
C 16:    p.    2g-33 
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c.  Frequency 

Frequency  may  be  measured  with  a  scanning 
superheterodyne  receiver,  but  this  method  has  a  low  POI  -for 
reasons  similar  to  those  o-f  using  a  rotating  antenna  -for  DF . 
A  channelized  receiver  provides  a  -frequency  measurement  but 
this  usually  is  a  much  more  expensive  system. 

Most  common  is  the  IFM  receiver  (see  section  on 
receivers  in  this  chapter),  which  has  an  excellent  resolution 
and  high  POI . 

d.  Pulse  Width 

Pulse  widths   can    be  measured    by    simply   noting 
when    the   pulse   rises    and   -falls    through    a    threshold. 
Re-flections    tend    to    severely    corrupt    this   method   o-f 
measur emen  t ,    however . 

A  more    accurate  method    is    to    sample    the    pulse  many 
times   and    take    the   PW    in    the    distance   between    the    3    dB 
points. 

e.  Time   o-f    Arrival     <TOA) 

Threshold   crossing  measurement    o-f    this    parameter 
could   also   be   used    as    in    the    case    o-f    the    PW,    but    a   better 
measurement    can    be   obtained   by   noting    the  TOA    o-f    the   -first    3 
dB   point.      TOA's   are   used    in    the   calculation    o-f    PRI's. 
Adjacent    pulse  TOA's   are    simply    subtracted    to    -find    the   PRI  . 
■f .       Higher    Order    Parameters 

Parameters    that    are   calculated   -from    the   other 
initial    parameters   are    the   derived    or    higher    order    parameters. 
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The  PRI  is  used  to  link  interrupted  or  split 
chains  o-f  pulses  -from  the  same  emitter.   It  can  also  be  used 
as  a  parameter  to  identify  a  jittered  or  staggered  PRI. 

Scan  patterns  can  be  recognized  by  analysis  o-f  the 
pulse  amplitudes  generated  as  the  radar  beam  sweeps  across  the 
platform.   Knowledge  o-f  the  scan  pattern  can  be  of  help  in 
identifying  the  emitter  since  the  radiation  pattern  can  be 
correlated  to  emitter  type.   [21:  p. 1673 
2.   Deinter 1 eaving 

A  stream  of  pulses  from  several  different  emitters 
includes  different  frequencies,  amplitudes,  and  PRI's  mixed 
together  (interleaved).   Separating  these  different  pulses 
into  chains  of  similar  pulses  for  identification  is  the 
deinter 1 eaving  problem. 

Ni th  high  density  environments,  the  time  available  for 
processing  each  pulse  may  be  only  around  1  microsecond.   C21: 
p. 168]   Therefore,  deinter 1 eaving  algorithms,  if  done 
digitally,  must  be  fast.   Analog  equipment  that  sort  pulses 
have  the  same  requirement.   This  required  speed  is  the 
obstacle  that  hinders  the  use  of   Digital  FFT's  for  pulse 
sorting. 

Provision  must  be  made  in  algorithms  for  missing 
pulses  in  the  data  stream.   Missing  pulses  is  one  of  the 
greatest  confusion  factors  in  the  system.   Methods  are 
currently  in  use  to  smooth  fluctuations  of  the  number  of 
pulses  that  are  input  into  the  preprocessor.   C21:  p. 1673 
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General  discussions  of  deinterl eaving  methods  are 
given  in  this  section.   The  -first  two  are  digital  methods, 
using  the  measured  parameter  of  the  pulse,  and  use  of  one  or 
the  other  depends  upon  the  accuracy  of  the  ESM  system.   Two 
hardware  methods  are  briefly  covered  also. 

a.  Cell  or  "Pigeon  Hole"  Techniques 

A  system  that  accurately  measures  the  parameters 
of  the  pulse  stream  can  use  this  technique  as  a  fast  and 
efficient  deinterl eaving  scheme. 

Pulse  data  with  similar  parameters  are  directed 
into  cells  or  "pigeon  holes".   The  cells  soon  contain  alike 
parameters,  such  as  bearing,  frequency,  or  pulse  widths  that 
can  be  analyzed  for  emitter  identification  or  further  sorting 
if  too  wide  of  a  variation  exists. 

Variations  could  be  caused  by  emitter  variations 
itself  or  purposeful  variations  such  as  a  jittered,  staggered 
PRI ,  or  frequency  agile  radar.   If  variations  (unmatched 
parameters)  are  present,  then  they  are  passed  to  a  more 
complex  algorithm  for  processing  for  the  special  modulations. 
Otherwise,  a  simple  analysis  or  comparison  with  the  active 
emitter  file  or  emitter  parameter  file  results  in  an  emitter 
identification . 

b.  Time  of  Arrival  Techniques  C21:  p. 168] 

When  accurate  computation  of  pulse  parameters 
isn't  available,  a  time  slice  of  the  pulse  stream  is  used  to 
simplify  the  data  under  analysis.   A  time  slice  will 
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contain  about  8  to  512  pulses.   This  scheme  reduces  the  number 
o-f  pulses  being  considered  in  a  high  density  enviroment,  and 
forces  comparison  o-f  parameters  to  those  only  in  the  present 
time  si  ice. 

The  -first  pulse  in  the  time  slice  is  used  as  a 
re-ference  and  the  rest  o-f  the  pulses  is  scanned  -for  a  match. 
When  a  match  is  -found,  the  TOA  o-f  the  pulses  is  used  to 
calculate  a  PRI  .   When  a  su-f-ficient  number  (usually  6   or  7)  o-f 
matching  pulses  are  -found,  then  the  data  is  grouped  as  a  chain 
and  is  output  to  the  main  processor.   A  problem  generated  with 
this  technique  is  the  splitting  o-f  chains  in  the  pulse  data 
stream  by  the  time  slice.   Several  time  slices  may  have  to  be 
analyzed  be-fore  enough  matches  in  a  slice  are  -found  to  -form  a 
chain . 

c.   Analog  Methods 

A  pulse  sorter  operating  on  the  pulse  period 
that  has  some  analogy  to  digital  TOA  methods  is  currently  in 
use.   C16:  2g-7  to  2g-9]    Consider  Figure  2.18. 

The  pulse  stream  is  -fed  into  a  coincidence 
detector-discriminator  that  attempts  to  match  the  pulses  with 
gates  generated  by  a  gating  pulse  generator.   The  gating 
pulses  are  initially  generated  at  a  PRI  slightly  higher  than 
the  highest  pulse  PRI.   By  gradually  reducing  the  gate  PRI, 
the  highest  PRI  pulse  train  will  be  captured  by  the  gates  i-f 
the  phase  di-f-ference  between  the  gates  and  pulses  o-f  the 
stream.   This  phase  di-f-ference  is  used  to  generate  an  error 
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signal  that  causes  the  gate  PRI  to  decrease  until  synchronized 
with  the  highest  PRI  pulse  train  o-f  the  stream.  The  pulse 
train  is  now  sorted  out  o-f  the  stream,  and  the  remaining 
pulses  are  applied  to  a  second  coincidence 
detector-discriminator  where  the  process  is  repeated. 

Other  methods  are  used  on  the  pulses  left  in 
the  stream  after  passing  throught  all  the 

detector-discriminators.   Staggered  PRI's  can  also  be  handled 
by  using  two  or  more  gate  generators  operating  in 
synchronization . 

d.   Frequency  Domain  Methods 

Bandpass   comb   -filters    C16:    p.2g-9]    can    be 
used    to    pass    only    harmonically    related   -frequency    components. 
By    tuning    the   -filter   with    the    pulse    stream    as    an    input,     the 
■filter   will    pass   components  when    it    is    tuned    to    a   PRI    o-f    the 
pulse   stream.       Disadvantages   o-f    this  method    include    the    low 
resolution    o-f    the    sorted   PRI's   and    inability    to   distinguish 
harmonic    components   o-f    different    PRI's. 

Digital    FFT's   are    attractive    as    a   sorting 
tool,    but    hardware    (size    and   cost)     and   computational    speed 
have   a  way    to   go   before    it    offers    practical    alternatives. 
VHSIC  will    certainly    improve   both    the   hardware    and   speed 
aspects. 
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Ill .   THEORY  OF  HALSH  AND  RADEMACHER  FUNCTIONS 

Electrical  Engineering  could  be  said  to  revolve  around  the 
sine  and  cosine  -functions.   They  are  the  basis  for  development 
in  many  areas  due  in  part  to  the  inherent  properties  in  their 
•frequency  domain  representations. 

The  application  of  digital  techniques  and  semiconductor 
technology,  though,  has  brought  -forth  uses  and  awareness  o-f 
other  orthogonal  functions.   These  o-ften  do  not  have  the 
desirable  properties  o-f  sine  and  cosine  -functions  -for  use  in 
linear,  time  invariant  networks,  but  they  do  have  other 
advantages  that  render  them  use-ful  in  other  applications. 

The  early  part  o-f  this  century  saw  the  introduction  o-f 
several  two  level  orthogonal  -functions,  or  sometimes  called 
binary  -functions  because  o-f  their  amplitudes  taking  only  two 
values.   Work  in  orthogonal  matrices,  by  Sylvester  in  1367 
and  Hadamard  in  1893,  was  an  early  approach  to  these 
■functions. 

In  1918,  Al -f red  Haar  presented  a  complete  set  o-f 
rectangular  -functions  that  took  only  two  values  but  yet 
provided  complete  expansion  o-f  a  continuous  function.   These 
Haar  functions  could  be  made  to  converge  uniformly  and 
rapidly.   See  Reference  23  for  additional  information. 

A  German  mathematician,  H.  Rademacher ,  presented  another 
set  of  orthogonal  functions  in  1922.   These  were  followed  by 
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the  Walsh  Functions,  de-fined  in  1923  by  the  American 
mathematician  J.L.  Walsh.   The  Rademacher  -functions,  although 
independently  presented,  were  -found  to  be  an  incomplete  but 
true  subset  to  the  Walsh  -functions.   C23:  p.v-vi] 

This  chapter  will  brie-fly  examine  the  Rademacher 
functions.   They  are  important  because  products  o-f  the 
■functions  yield  a  certain  ordering  o-f  the  Walsh  -functions. 
The  Walsh  functions  will  then  be  examined  in  more  detail. 
The  reader  is  referred  to  Appendix  B  for  a  generic  treatment 
of  orthogonality  and  orthogonal  functions. 

A.   RADEMACHER  FUNCTIONS 

These    two-level    orthogonal     functions   are    represented    by 

R<n,t>  <3.1) 

n    =    8,1,2,3,... 
and   can    be   seen    in    figure    3.1    for    n   =   8    to    6.       [23:    p. 61 

Notice  that  they  are  a  series  of  rectangular  pulses  or 
square  waves  and  have  2n_1  periods  over  a  time  base  from  8 
to  T.      Mathematically    they   could   be   defined   by 

R<n,t>    =    sign<sin<2n*t>>  <3.2) 

The    first    Rademacher    function,    R<8,t>,    is   equal     to    one    for 
the   entire    interval ,    and   subsequent    Rademacher    functions   have 
odd   symmetry  with    amplitudes   of    +1    and    -1.      All    have   unit 
mark-space   ratio. 

One   can    generate    these    functions  with    a    sine    function    of 
appropriate   frequency  with    amplification    and   hard    limiting. 
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The  appropriate  frequency  would  be  one  with  the  same  zero 
crossing  positions  as  the  Rademacher  functions.   [23:  p. 73 


R(0,D 


-I 


R(l.t) 


R(2,t) 


R(3,l) 


.R(4,t) 


m 


.R(5.l) 


Figure   3.1.      The   First    6   Rademacher    Functions. 

B.      WALSH    FUNCTIONS 

Of   more    importance    than    the   other    two-level    orthogonal 
functions   discussed    are    the  Walsh    functions.      They   form    a    set 
of    rectangular   waveforms    taking   only    two    amplitudes,    +1    and 
-1,    but    do   not    have   a   unit   mark-space   ratio    like    the 
Rademacher    functions.      They   are   defined   over    a    time    interval 
T,    which   must    be   known    to    assign   values    to    the   functions. 

Two   arguments   are   required   for    complete    specification    of 
a  Wal sh    function 

WAL<n,t>  (3.3) 
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The  time  base  Ct]  of  the  -function  is  usually  specified  as 
t/T,  and  thus  normalized  -from  9  to  1.   The  number  Cn]  is 
equal  to  the  number  o-f  zero  crossings  a  Walsh  function  has 
during  the  time  base.   C23:  p. 73 

Figure  3.2  shows  the  Walsh  functions.   The  symmetry  of 
the  functions  and  the  concept  of  sequency  (next  section)  was 
used  by  Harmuth  C243  to  define  another  notation  for  each 
Walsh<n,t).   Each  WAL(n,t)  is  either  odd  or  even  about  the 
midpoint.   If  the  function  is  odd,  then  it  can  be  referred  to 
as  SAL<k,t).   If  the  function  is  even,  then  WAL<n,t)  is  also  a 
CAL<k,t).   Note  the  similarities  between  this  notation  and  the 
sine-cosine  functions.   Relationships  between  Cn]  and  Ck]  are 

CAL<k,t)  =  WAL(2k,t)  <3.4) 

SAL(k,t)  =  WAL<2k-l,t) 
This  notation,  SAL<k,t)  and  CAL(k,t),  is  seen  on  the  right 
side  of  Figure  3.2. 

1.   Sequency,  Ordering,  and  Phasing 

Sequency  is  a  term  proposed  by  Harmuth  to  describe  a 
periodic  repetition  rate  which  is  independent  of  waveform.   It 
is  defined  as  "one  half  of  the  average  number  of  zero 
crossings  per  unit  time  interval"  and  abbreviated  "Zps" .  C23: 
p.  133   A  188  HZ  sine  wave  has  288  zero  crossings  per  second, 
so  frequency  is  a  special  measure  of  sequency  as  applied  to  a 
sinusoidal  waveform. 

Ordering  the  functions  by  ascending  number  of  zero 
crossings  <WAL<n,t)  notation)  is  called  sequency  ordering. 
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Figure  3.2.   The  Walsh  Functions  in  Sequency  Order 
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Although  Cn]  is  not  the  sequency  of  the  -function,  the 
equivalence  o-f  the  WAL(n,t)  with  SAL(k,t>  and  CAL(k,t) 
notations  automatically  arranges  the  functions  in  ascending 
order  o-f  sequency.   The  value  o-f  Ck3  in  the  SAL  and  CAL 
-function  is  the  sequency  -for  the  -function.   Sequency  ordering 
is  the  preferred  ordering  -for  spectral  analysis  and  -filtering 

The  natural  ordering  (or  Paley  ordering)  is  obtained 
by  generating  the  Walsh  -functions  with  products  o-f  Rademacher 
■functions.   Figure  3.3  shows  the  Walsh  -functions  in  natural 
order,  which  are  referred  to  as 

PAL<n,t)  <3.5) 

n    =    @,1,2,3,... 
because   Paley    [24]    used    this    ordering.       For    theoretical     and 
mathematical    work,    image    transmission,    and    computational 
efficiency,    this   ordering    is   usually    the   preferred   one.       C23: 
p.     18] 

A  modulo-2    addition    (Appendix    C)     relationship    exists 
between    the  Walsh    functions  with    sequency   ordering   and    the 
functions  with    natural    ordering.       C23:    pp. 31-32]      An    example 
will     illustrate    the    conversion    between    the    two    order ings. 
Consider   WAL(nw,t>    =  WAL(13,t>. 

Here    nw  =    13  =    1 18 1 2 -      To    find    the    equivalent 
natural    order    function,    PAL(np,t),    first    add   a   zero    (base 
2)    as    the    leftmost    digit    to    base   2   nw.      Then    generate    the 
value   of    np    in    base   2   by   doing  modulo-2    addition    between 
consecutive   zeroes    and    ones   of    base   2    nw,    starting  with    the 
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Figure   3.3.      The  Walsh    Functions    in   Natural    Order 
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leftmost  digit  (the  zero  that  was  added  on  to  base  2  nw)  . 
Following  these  rules,  nw  =  13  =  H8I2  =  911912?  and 

8  $  1  =  1  =  k3 

1  e  1  =  8  =  k2 
1  $  9  =  1  =  kt 

9  $  1  -  1  =  k0 
So, 

np  =  <k3k2k1k9)2  =  19112  =  11 

Figure  3.3  shows  the  result  is  correct. 
WAL( 13,  t)  =  PALC  11, t) 

The  phasing  o-f  the  -functions  as  given  in  Figure  3.2 
and  3.3  is  known  as  Harmuth  [25]  phasing.  This  phasing 
emphasizes  the  phase  similarities  with  sine-cosine  functions. 

Positive  phasing  [23:  p. 18]  is  where  all  of  the 
functions  begin  at  +1.   This  involves  a  sign  change  for  many 
of  the  Figure  3.2  and  3.3  functions,  and  is  the  result  of 
function  derivation  from  Hadamard  matrices  (orthogonal 
matrices  composed  only  of  + 1  and  -1  elements;  see  Ref .  23, 
pages  24-25)  or  Rademacher  functions  (next  section). 
2.   Derivation  of  the  Walsh  Functions 

Products  of  selected  Rademacher  functions  will  yield  a 

complete  set  of  Walsh  functions  in  natural  order.   Recall  that 

the  natural  ordering  of  the  Walsh  functions  was  referred  to  as 

PAL(n,t) . 

In  terms  of  Rademacher  functions, 

m+  1 
PAL(n,t)  =r~[~P  btR(i  ,t)  (3.6) 

i=l 

where  b^  is  either  a  zero  or  a  one,  indicating  the  presence 
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or  not  o-f  the  i'th  Rademacher  -function.   To  find  the  bj's 
and  Cm],  represent  Cn]  as  a  binary  number, 

»  =  <bm+lbmbm-r--b1)2         <3-?> 
with  m  =  the  highest  power  o-f  two  -found  in  the  binary  number 

Cn].   The  bj's  present  in  the  positions  o-f  base  2  Cn] 

indicate  the  presence  or  not  o-f  the  Rademacher  -function  in 

Equation  (3.6) . 

Consider  PAL(?,t) .   Then 

n  =  9 

n  =  18912 

with  m  =  3,  b4  =  1,  bj  =  1,  and  -finally, 

PAL(?,t)  -  R<4,t)R<  l,t) 

It  should  be  noted  that  this  product  is  positive 
phasing.   Harmuth  phasing  can  be  obtained  by  de-fining  the 
Rademacher  -functions  over  the  interval  -34<t<^.   This 
de-finition  inverts  all  o-f  the  -functions  shown  in  Figure  3.1. 
Products  o-f  the  inverted  -functions  will  yield  natural  ordered 
and  Harmuth  phased  Walsh  -functions.   C23:  p. 22]   The  Walsh 
•functions  with  sequency  ordering  can  be  derived  with 
Rademacher  -functions  by  use  o-f  the  Gray  Code  (Appendix  C)  . 
Consider  WAL(13,t)  where  n  =  13. 

In  Gray  Code,  n  =  13  is  represented  by  1811,  and 
recognizing  the  ones  of  the  Gray  Code  are  in  the  fourth, 
second,  and  first  positions,  (similar  to  the  derivation 
of  the  PAL(n,t)  functions),  C23:  p. 23] 

WAL(13,t)  -  R(4,t)R(2,t)R(  l,t) 
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Other    derivations   of    the  Walsh    functions    can    be 
accomplished   by   di-f-ference    equations,    Hadamard  matrices,    and 
Boolean    synthesis.       See    Re-ference   23,    pages   28-26. 
3.      Walsh    Series 

It    is  well    known    that    a    time   -function    -f(t)     can    be 
expressed   as    a   sum   o-f    a   series    o-f    sine   and   cosine   -functions. 
Each    -function    has    a    coe-f -f icien t    that    determines    the  value    o-f 
the   -function    in    that    series. 

The    time    -function    can    be   expressed    in    a   similar    way 

using    the  Walsh    -functions.       C23:    p.  13] 

N-l 
f<t>    =    a9WAL<8,t>     +      Z      anWAL<n,t)  <3.3) 

n=l 

where  the  coefficients  would  be  calculated  from 


J8 


ae/2  =  1/T  ■  \  f  <t)WAL<8,t)dt  (3.9) 


J8 


an  =  1/T  •  \f  <t)WAL(n,  t)dt  (3.18) 

J8 

with  N  =  the  desired  number  of  Walsh  terms,  and  T  is  a  given 

length  of  the  function. 

Each  Walsh  function  has  a  coefficient  associated  with 
it  that  gives  the  value  of  the  function  in  the  series  the  same 
way  as  does  the  Fourier  Series  coefficient.   C23:  p.  13] 

If  the  function  is  periodic,  and  T  is  normalized  to  be 
equal  to  1,  then  a^/2  is  the  mean  of  the  function,  since 
WAL<8,t>  is  equal  to  1  over  the  interval  8<t<l.   The 
expressions  for  ag/2  in  both  the  Walsh  and  Fourier  Series 
are  equal  with  this  normalization. 
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Consider  that  the  -function  is  absolutely  integrable 
in  the  interval  8  < t  <  1 .   Then  the  -function  can  be  expanded  in  a 
Nalsh  series  of  the  -form   C26:  p. 232  and  23:  p. 48] 

f<t)    =    Z   a<n)WAL<n,t)  <3.11) 

n=8 


where 


"■-$: 


a<n)    =    \    f  <  t)WAL<n,t)dt  <3.12) 

Je 

A  more    general     definition    -for    a   periodic    -function    defined   over 
an    interval     (8,T)     is 

04 

f<t)    -    Z   a<n)WAL<n,t/T)  <3.ila) 

n=9 


a<n)    =    1/T    )    f  <t)WAL<n,t/T)dt  <3.12a) 

8 

Now   coefficient    ag    becomes    the   mean    of    the    function, 
since  WAL<8,t/T)    =    1    over    the   normalized    interval     8<t/T<l. 

The   Walsh    Series   can    be   defined   over    the    same    interval 

<8,T>    using    the   CAL    and   SAL   notations,    and    an    infinite    number 

of    terms:       E26:    p. 2323 

oo  oo 

f<t)    =    a0WAL(8,t)     +       Z  Z    < a^ SALC i , t)     +    bjCAL<j,t)) 

i=l         j=l 

<3.  13) 

wi  th 

as<n)        1/T        \f  <t)SAL<n,t/T)dt  3.14) 


9. 


ac<n)    =    1/T    •    \  f  <t)CAL<n,t/T)dt  <3.15) 
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A    couple    of    simple   examples   -follows    to    illustrate    the 
Walsh    Series   expansion. 

a.       Expansion    o-f    Sin(t)     in    a  Walsh    Series 

The    sine    -function   will    be    de-fined   over    one   period, 

with    a   period    o-f    2*,    and    the  Nal  sh    functions  will    be    de-fined 

as    t/T   goes   -from   8    to    1 . 

Calculating    the    -first    coe-f  f  icient ,    using   Equation 
<3. 12a) , 

p2* 
a9    -    1/2*    ■    \sin<t)WAL(9,t/T)dt 

p2* 
a0   =    1/2*    •    \sin<t)  <  1)  dt 

a0   =   8 

This    is    a    con-fident    calculation,    as  most    students 

Know    that    the    average    of    a    sinusoid    is    8! 

The   second    coefficient, 

p2* 
at   -    1/2*    •    \  sin<t)WAL<t,t/T)dt 
J  9 


•    C \sin 
J9 


n*  p2* 

at   =    1/2*    •    C\sin<t)  <-l)dt    +      \sin<t)<l>dt    ] 


aj   -   -2/*     =      -.63611 

A  lot  of  menial  integration  will  yield  only  these 

coefficients  of  any  practical  value: 

a5  =  .26348 
a13  =  . 12653 

Coefficients  numbered  2,6,9,  and  18  are  present 

but  have  values  less  than  8.86. 
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The   -function    sin<t)     is    then    approximated   by 
sin<t>    =      a1l»iAL<lft)     +    a^MAL<5,t)     +    a  13WAL< 13 , t> 
where    the   coefficient   values    are    given    above. 

Figure    3.4    C23:    p.     153    shows    the    three  Walsh 
■functions    used    and    their    sum.      The   negative    sign    o-f    the 
coefficient    a^    has   been    used    to    invert    the  Walsh<l,t) 
function . 


,WAL(I,I) 


.WAL(5,t) 


ir~ i_nj~L^ 


r& 


WAL(l3,t) 


ZWAL(l,t),WAL(5,ll,WAL(l3,t) 


Figure  3.4.   Sin<t)  with  a  Limited  Nalsh  Series 
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Of  course,  the  addition  o-f  more  of  the  terms  of 
the  Series  would  yield  a  better  approximation.   One  can  see 
that  representing  a  smooth  curve  using  rectangular  shapes 
would  require  a  large  number  o-f  terms  o-f  the  Walsh  Series. 

The  Walsh  Series  will  give  better  results 
representing  a  rectangular  function. 

b.   Expansion  of  a  Rectangular  Function 

A  rectangular  pulse  function  should  lend  itself 
nicely  to  representation  by  other  rectangular  functions. 

Indeed,  it  is  hoped  that  the  representation  might 
be  exact  or  with  small  error  in  only  a  few  terms. 

Consider  this  function: 

f<t)  =  1  8<t<.25 
8  .25<t<.75 
1      .75<t<1.8 

Proceeding  with  the  coefficients, 
a0  =  \f  <  t)WAL<8,t)dt 


■i 


:  25         ri .  8 

(D(l)dt   +  \   (D(l)dt 
8  J. 75 


and, 


=  8.5 


\ 25  r  1 . 8 

<  1)  <-l)dt 


-  \  <  1)  (-l)dt   + 


75 


=  -.5 
Coefficients  a^,  ag,  and  above  will  all 
calculate  to  be  zero. 
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This  rectangular  time  -function,  which  could 
represent  two  radar  pulses,  can  be  represented  exactly  by  two 
Walsh  functions. 

f<t>  -  8.5WAL<8,t>  -  8.5WAL(2,t> 

Figure  3.5  shows  the  -function,  the  Walsh 
-functions,  and  the  sum  o-f  the  Walsh  -functions.  This  is  quite  a 
di-f-ference  -from  the  representation  by  Fourier  Series,  which 
would  require  a  lot  o-f  terms  to  present  an  inexact 
representation . 

c.   Wave-form  Synthesis 

These  results  obtained  above  indicate  that  a 
continuous  wave-form  is  more  suited  -for  Fourier  trans-formation 
and  a  discontinuous  wave-form,  certainly  a  rectangular  one,  is 
more  suited  for  Walsh  transformation.   This  conclusion  should 
be  easy  to  accept  with  examination  of  the  structure  of  the 
orthogonal  functions  (sine-cosine  and  Walsh). 

Figure  3.6  C23:  p.  333  shows  a  reconstruction  of 
two  waveforms  and  the  number  of  Fourier  and  Walsh  terms  used. 
The  reconstruction  supports  the  conclusion  ^ery   well.  In 
addition,  Beauchamp  C23:  p. 353  reports  that  the  same 
conclusion  can  be  drawn  when  the  Fourier  and  Walsh  Transform 
are  used  to  reconstruct  a  continuous  and  a  step  waveform. 

In  coming  to  this  conclusion,  Beauchamp 
considered  the  effect  of  the  number  of  allowed  levels  for 
quantizing  the  waveform  and  the  sampling  interval.  C23:  p. 35] 
The  sampling  interval  determines  the  number  of  coefficients 
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L 

Figure  3.5.   Expansion  o-f  a  Rectangular  Function 
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Origmol 


Walsh    131  terms 


4   samples  puis* 


A 


Origmol 


Wolsh  2  4   terms 


Fourier    18  terms 


Fourier    44   terms 


*-«>      \l  >* 


Figure    3.6.      Wave-form   Synthesis   Using  Walsh    Functions 
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produced,  and  the  highest  Walsh  coe-f -f icient  -found  -from  the 
Trans-form  determines  the  number  o-f  data  points  available  -for 
any  quantized  level.   This  is  a  built-in  limitation  on  the 
accuracy  of  the  Walsh  Trans-form  that  the  Fourier  Trans-form 
doesn't  have. 

d.   Digital  Sampling 

A  theorem  similar  to  the  Sampling  Theorem  -for 
sine-cosine  functions  is  applicable  to  the  Walsh  -functions. 
The  minimum  sampling  rate  is  C23:  p. 373 

■fs  =  2k+1  <3.  16) 

with  k  being  the  power  o-f  2  that  represents  the  sequency 
bandwidth.   The  sequency  bandwidth  would  be  the  sequency  o-f 
the  highest  component  desired  to  be  represented. 

Suppose  a  CAL(12,t)  component  had  been  determined 
and  it  was  desired  to  include  this  component  in  the  bandwidth. 
The  sequency  o-f  this  component  is  12,  and  the  closest  power  o-f 
2  that  is  at  least  12  is  4  <24  =  16)  .   Thus  k  =  4  and 

fs  =  24+1  =  32 
For  comparison,  a  cosine  o-f  -frequency  12  would 
require  a  sampling  rate  o-f  24  <2*-frequency)  .  Aliasing  is 
equally  applicable  to  the  Walsh  -functions  and  should  be 
considered.   £23:  p. 373 

4.   The  Walsh  Trans-form 

Section  3  stated  that  a  continuous  -function  over  the 
interval  8<t<l  could  be  represented  by  a  sum  o-f  Walsh 
•functions,  shown  in  Equation  3.11. 
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M 

-f<t)    =      Z   anWAL<n,t)  <3.11) 

n=8 


where 


-5 


<t)WAL<n,t)dt  (3.12) 


Thus,    a    trans-form   pair    can    be    de-fined:    C23:    p.  48] 

M 

-fCt)    =      Z   F(k)WAL(k,t)  <3.17) 

n=8 


K)    =\ 


F<k)    =  \    -f<t)HAL<k,t)dt  <3.13> 

'8 

F<k)  could  be  written,  -for  an  interval  8<t<T, 

F<k>  =  1/T  \  -f(t)HAL<k,t)dt        (3.1?) 


j: 


Now  suppose  the  interval  T  is  divided  into 
N  parts,  letting  t  =  iJi,  where  -J  is  a  small  time  increment  and 
i  is  an  integer.  The  aim  is  to  convert  the  integral  into  a 
summation  -for  discrete  computation. 

Then  , 


Je 


wi  th 


F<k>  =  1/T  \  -f<  t)WAL<k,t)dt 


N-l 
%    l/NA    •  Z  -f<i^l)WAL<k,i4> 
i=8 


A   ■  dt 
NA   =  T 
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Remembering    A    to    be    a   given    small     increment,    but 

unchanging    in   value, 

N-l 
F<k)    -    1/N    •    Z   f<i)WAL(k,i)  (3.22) 

i=8 

wi  th 

f<i)    =   evaluation    o-f    f<t)    at    intervals    i  A 

=      a   series   o-f    numbers 

Change    the    notation    from   k    to   n,    -f<i)     to 

x  j  ,    and   F<k)    =  Xn ,    and  write    the   -finite   discrete  Walsh 

trans-form   pair. 

N-l 
Xn   =    1/N    ■    Z  XiWAL<n,i)  < 3 . 23) 

n=9 

n   =    8,1,2,.. .N-l 

and, 

N-l 
Xi    =         Z  XnWAL(n,i)  < 3 . 24) 

n=8 

i    =    8, 1,2, .. .N-l 

The  CAL  and  SAL  trans-forms  for  a  discrete  series  can 

also  be  found  by  using  the  equivalence  of  WAL  with  CAL  and  SAL 

notations: 

N-l 
Xc(k)    =    1/N    •    Z  XiCAL<k,i)  <3.25) 

i=8 

N-l 
Xc<k)    =    1/N    -    Z  XjSAL<n,i)  <3.26) 

i=8 

An    even    or    odd    sequence  when    transformed  will    have 

similar    properties   as    the   discrete   Fourier    Coefficients    from 
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even  or  odd  functions.   That  is,  an  even  sequence  (symmetrical 
around  its  midpoint)  will  trans-form  into  only  CAL  -function 
coe-f-f icients,  and  an  odd  sequence  will  trans-form  into  only  SAL 
coe-f-f icients.   C23:  p.41] 

The  Discrete  Walsh  Trans-form  <DWT)  has  some  advantages 
over  the  Discrete  Fourier  Transform  (DFT)  because  it  involves 
only  additions  and  subtractions,  not  multiplications  as  in  the 
DFT.   Also,  the  transform  is  not  noisy  because  precise 
representations  in  the  digital  computer  is  possible.   Precise 
representations  of  sine-cosine  functions  in  the  DFT  is  not. 
Finally,  the  Walsh  transform  kernel  is  ±1  and  is  considerably 
easier  to  calculate.   C23:  p. 423 

An  important  point  of  the  transform  is  that  the  DWT  is 
its  own  inverse,  and  a  separate  inverse  transform  definition 
is  not  required.   See  Equations  (3.23)  and  (3.24). 

a.   The  Fast  Walsh  Transform 

Computation  of  the  DWT  would  involve  N* 
mathematical  operations  of  either  addition  or  subtraction. 
This  number  of  operations  has  been  decreased  by  using  the 
redundancy  in  the  Walsh  matrix  representation  of  the  Transform 
to  Nl og^J  operations,  the  same  number  as  in  the  DFT.   C23: 
p. 543   Remember  that  the  DWT  involves  additions/subtractions, 
not  complex  multiplications/additions  as  in  the  DFT. 

The  Fast  Walsh  Transform  can  be  described  in  a 
signal  flow  diagram  that  results  in  a  "butterfly"  appearance 
just  as  in  the  FFT.   See  Figure  3.7  on  the  next  page.    Unique 


88 


Figure   3.7.      FWT    "Butter-fly":    N  =    16 
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■features  of  the  butter-fly  in  the  -final  stages  o-f  calculations 
o-f-fer  an  opportunity  to  calculate  the  FWT  "in-place",  without 
using  computer  working  space.  C23:  p. 55]   The  result  is  not  in 
sequency  order,  however,  and  a  sorting  routine  must  be 
included  i-f  sequency  order  is  desired. 

Appendix  A  gives  listings  o-f  the  FWT  programs  used 
in  this  thesis. 

b.   Walsh  Coe-f-f  i  cients  and  the  Walsh  Matrix 

The  computation  o-f  the  coe-f-f i cients  -from  the  FWT 
is  \/ery    fast  and  takes  advantage  of  the  properties  of  a 
Hadamard  matrix   C23:  p. 52] 

The  Walsh  matrix  is  related  to  the  Hadamard  matrix 
and  can  be  used  to  determine  the  coefficients  of  the  Walsh 
Series.   C27:  p. 4]   Consider  this  matrix  equation, 

Can]  =  1/N  •  EW]  •  lxil  <3.27> 

where  CW]  is  a  N  X  N  matrix  of  the  Walsh  functions,  Cxj]  is 
a  column  vector  of  N  sampled  values  of  the  function,  and 
Can] ,  n  =  9  to  N- 1 ,  is  a  column  vector  of  the  Walsh 
coefficients. 

The  Walsh  matrix  can  be  thought  of  as  the  matrix 
of  sampled  Walsh  functions,  and  thus,  is  composed  of  plus  and 
minus  ones.   To  compute  the  coefficients,  simply  perform 
addition  and  subtraction  of  the  sampled  values  according  to 
the  row  of  the  Walsh  matrix  being  used.   An  example  will  show 
this  computation. 
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Let  the  number  o-f  samples  in  the  column  vector  =  3 
=  N.   Nriting  the  matrix  equation,  omitting  the  1/N  term, 


Can3  =   CW] 


Cxt3 


(3. 16) 


a8 

M 
a2 

*3 
a4 

a5 
a6 
a7 


+  +  +  +  +  +  +  + 

+  +  +  +  ---- 

+  ♦-  —  -  —  +  + 

+  +  --  +  +  -- 

+  --  +  +  --  + 

+  --+-  +  +  - 

+  -  +  --+-  + 

+  -  +  -  +  -  +  - 


x8 

xl 

x2 

x3 

x4 

x5 

x6 

x7 

For  comparison,  use  N  sampled  values  o-f  the 
rectangular  -function  on  page  73.   Thus, 

Cxi3T  =[11988911] 

i   =   8,1,2,3,...,' 

Now   evaluate    the   constants   an    using  matrix 
multiplication.      A   plus   sign    in    the  Nal sh   matrix    equals   +1, 
and   a  minus   sign   means   -1.       Don't    -forget    to    divide   each 
constant    by  N. 


ag   =   8.5 
at   =   8 
a2  —   8-5 
a3  =   8 


a4  =  8 
a5  =  8 
a6  =  8 
ay  =   8 


Comparing  with  the  value  o-f  a2  given  in  the 
Nalsh  Series  expansion  o-f  this  same  -function  (page  73),  we 

note  a  sign  change.   The  reason  lies  in  the  use  o-f  positive 

phasing  -for  the  Walsh  matrix  and  the  use  o-f  Harmuth  phasing  in 
the  Series  expansion  computation. 

The  same  answer  has  been  reached,  although  this  is 
the  discrete  case. 
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Matrices   can    also   be   used    to    show    the    Inverse 
Discrete  Nalsh   Trans-form.       From    Equation    (3.24), 

N-Cx^    -    CW3 -Can]  <3.28> 

=    CW3 -CW3 -Cx^ 
and 

Cxt3    =    1/N-CW3 -CW] -[Xj]  <3.2?) 

Equation    (3.2?)     says    to    get    the    sampled  values   o-f 
the   original    within    a   constant   N,    just   multiply    the 
coe-f-f icien ts   by    the  Walsh    matrix    CW]  . 

c.       Input   Time   Shi-fts    and  Their    Effect 

The   DFT  magnitude    is    invariant    to    the   phase   o-f    the 
input    signal.      Although    Walsh   Trans-form    signals    conform    to 
Parseval's  Theorem,     the    DWT    is   not    invariant    and    the    same 
spectral    representation    cannot    be   achieved    independently    o-f 
the   phase   or    time   shi-ft    o-f    the    signal.       C23:    p. 423 

The   DWT    is    time    invariant   when    the    time    shift    is 
obtained   by    dyadic    translation    C23:    p. 433.      Namely,    if    Xj, 
the    input    number    series,    is    shifted    to   Zj    =    x,,^,    where 

i$p  =  modulo-2    addition 
then  , 

ZC2<K>    +    Zs2<k)    =  XC2<K)    +   Xs2(k) 
if    the   coefficients    are    expressed    in    CAL   notation    <ZC,    Xc) 
and   SAL   notation    <ZS,    Xs) .       See  Appendix    C   for    modulo-2 
addi  tion . 

The    importance    of   variations  with    time   shift    can 
be    lessened   by   considering    sums   of    squares   of    coefficients. 
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This  is  one  definition  o-f  the  Walsh  Power  Spectrum,  and  is 
considered  in  the  next  section. 

5.   The  Walsh  Power  Spectrum 

Spectral  analysis  using  sequency  has  some 
advantages  over  analysis  in  the  frequency  (Fourier)  domain. 
A  de-fined  power  spectrum  -from  the  Walsh  -functions  is  easily 
applied  to  discontinuous  or  time  limited  -functions,  for 
example,  and  it  is  possible  to  have  a  sequency  limited 
spectrum  for  a  time  limited  function  in  contrast  to  the 
Fourier  Spectrum.   C23:  p. 87] 

Although  there  are  several  derivations  of  the 
Walsh  power  spectrum  C23:  p.? 13,  a  more  normal,  i.e. 
analogous  to  Fourier,  approach  will  be  used  to  make  a  spectrum 
definition.   As  there  are  readily  available  methods  (FyMT)  for 
calculating  the  Walsh  coefficients,  a  combination, 
specifically  the  sum  of  the  squares  of  related  coefficients, 
would  give  an  easy  to  use  definition.   Also,  as  stated 
earlier,  by  using  the  sum  of  the  squares  of  certain 
coefficients,  the  time  shift  variance  that  plagues  the  Walsh 
Transform  can  be  lessened. 

Define,  then,  the  Power  Spectral  Density  (PSD) 
coefficients  of  the  Walsh  spectrum.   [23:  p. 188] 

P(8)  =  Xc2(8,t)  (3.31) 

P(k)  =  XC2(K)  +  Xs2(k) 


P(N/2)  =  Xs2(N/2) 


k  =  1,2,... ,(N/2  -  1) 
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Note    that    (N/2   +    1)    spectral    points    are 
generated.      This  method    is   analogous    to   Fourier    power    spectral 
analysis,    where    the   power    coefficient    is    the   sum   of    the 
squares   o-f    the   real    and    imaginary   parts   of    the   complex    Fourier 
Trans-form   coeff icients. 

Beauchamp    C23:    p. 188]    states    that    taking    the 
square   root    o-f    POO    does   not    give   a   sequency    amplitude 
spectrum,    but    other    re-ferences    126:    p. 233]    use    -P^    as    the 
definition    o-f    the   sequency   spectrum. 

Another    spectrum   can    be   de-fined   using    the  Walsh 
Coe-f-f  icients.      The   name   Group   Spectrum    is    the   author's    term 
and  will    be   used   here.       Beauchamp    C23:    p.  186]    refers    to    it    as 
"the   odd-harmonic    sequency    spectrum",    and   Campanel 1  a/Robinson 
C26:    p. 234]    reier    to    this   spectrum   definition    as    a   generalized 
Walsh-Fourier    spectrum    [23]. 

The   power    content    of    individual    groups   of 
frequency   components   has   been    shown    to   be   equal     to   groups   of 
components    in    the   sequency   domain.      Any   group   of    components, 
in    either    domain,    is   a   collection    of    the   orthogonal    components 
that   make   up    the   original    signal    or    function. 

Call    groups   of    sequency   components   Gn ,    with    En] 
an    integer,    and   consider    a   discrete   signal    made   up   of    a 
sequence    <sjj>    of    length   N. 

The    total    energy   of    the   signal     is,    by   Parseval's 
Theorem,    equal     to    the   sum   of    the   orthogonal    components    C29]. 
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<These  orthogonal  sequency  components  have  been  grouped  into 

6n.)   Then, 

N-l       m 
1/N  •  Zsk2  =  2Gn  <3.32) 

k=e *    n=e 

and  the  number  o-f  groups  Cm]  is  given  by 

m  =  1  +  1 og2N  <  3 . 33) 

The  groups  Gg  and  G<log2N)  represent  the  power 
content  o-f  the  d.c.  component  and  the  -folding  frequency 
component,  respectively.   The  other  Groups  represent  the  power 
content  of  the  component  at  n  =  1 ,  n  =  2,  n  =  3,  etc.,  plus 
all  odd  harmonics. 

G8  =  P2<8)  <3.34) 

Gi    =    P2<  i)     +    P2<3)     +    P2<5)     +     ... 

G2   =    P2<2)     +    P2<6)     +    P2<  19)     +     ... 

G3   =    P2<4)     +    P2< 12)     +    P2<28)     +     ... 

G4   =    P2<8)     +    P2(24)     +    P2(48)     +     ... 
and   so   on,    until 

G<log2N)    -   P2<N/2) 
The   plot    of    Gn   versus   n    gives    a    discrete    spectum 
that    is   unique    for    a    certain    time    sequence    C26:    p. 234],    and 
invariant    to    a    time   shift    of    the    input    C23:    p. 196]. 

This   spectrum    is   highly   compressed    <m   points 
instead   of   N/2   +    1    points),    and    there-fore   doesn't    present    all 
characteristics   of    the    data    as  well    as    the    previously    defined 
spectrum.       C23:    p. 186] 
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6.   The  Application  of  the  Fast  Walsh  Trans-Form 
The  application  o-f  the  Walsh  Transform  to 
continuous  or  discontinuous  -functions  is  easily  accomplished 
using  the  FWT  programs  provided  in  Appendix  A.   The  programs 
input  a  number  series  (which  could  represent  a  sampled 
wave-form)  and  compute  the  FWT.   With  these  coe-f-f icien ts,  the 
Power  Spectral  Density  and  the  Group  spectrum  can  then  be 
calculated.   The  FWT  will  be  applied  to  the  previous  examples 
and  the  results  will  be  compared  and  analyzed  in  terms  o-f  the 
elaborated  theory  in  this  chapter. 

a.   The  FWT  o-f  a  Sinusoid. 

Figure  3.8  gives  a  normalized  graph  o-f  the 
coe-f-f  icien  ts  o-f  the  FWT  o-f  one  cycle  o-f  a  sinusoid  in  the 
Walsh  interval  8  to  1.   The  actual  magnitudes  can  be 
determined  -from  the  scale  listed  on  the  vertical  axis. 
Running  either  program  in  Appendix  A  will  give  a  listing  o-f 
the  coe-f-f  icien  t  magnitudes. 

Note  that  the  maximum  coe-f-f  icient  (normalized 
to  1)  is  a  WAL(n,t)  =  WAL(l,t>.   Equation  <3.4)  o-f  the 
chapter  reveals  there  is  a  SAL(l,t)  with  a  sequency  o-f  1  at 
this  Walsh  index.   Thus,  the  one  cycle  sin  wave  in  the  Walsh 
interval  8  to  1  has  a  distinguishing  feature  at  a  sequency  of 
1.   Perhaps  not  surprising  since  it  is  known  that  frequency  is 
a  special  measure  of  sequency,  and  the  one  cycle  in  this 
interval  can  be  thought  of  as  having  a  frequency  of  1. 


88 


CO 


<: 

CO 

oc 

h- 

LL. 

a 

ZC 

co 

LU 

.J 

_J 

< 

u 

2e 

>- 

cj 

\- 

co 

*— 1 

< 

LL. 

CD 

CD 

GO 

SI 

•<* 

C 

CM 
CO 

UJ 
Q 

*-H 

CM 

X 

CO 

-J 
< 

CD 


CO 


s 


in 


s 


in 


(1969  = I)    OVW  138 


Figure  3.8.   FWT  of  1  Cycle  of  a  Sine  Function 


8? 


Also  note  the  other  prominent  components  at 
n  =  5,  9,  13,  and  29,  with  sequencies  of  3,  5,  7,  and  15, 
respectively.   Our  previous  calculations  with  the  Walsh  series 
(page  71)  predicted  conspicuous  coeff icien ts  at  n  =  5  and 
13. 

Applying  equations  (3.31)  to  these 
coe-f-f  icien  ts  results  in  Figure  3.9,  the  Power  Spectral  Density 
(PSD)  .   The  maximum  value  of  power  is  at  a  sequency  of  1,  and 
there  is  no  DC  power  term  (sequency  of  9) .   Harmonic  power  is 
located  at  sequencies  of  3,  7,  and  15. 

These  odd  number  power  components  combine 
into  a  single  component,  Gj,  in  the  Group  spectrum,  Figure 
3.18.   See  Equation  (3.34). 

An  increase  of  frequency  gives  similar 
results.   Figures  3.11  to  3.13  show  the  FWT ,  PSD,  and  Group 
Spectrum  of  three  cycles  of  a  sinusoid  in  the  Walsh  interval  0 
to  1.   The  maximum  FWT  coefficient  lies  at  n  =  5  Ca  SAL(3,t) 
component  with  a  sequency  of  3],  with  other  major  components 
at  n  =  1,  9,  13,  21,  25,  29,  53,  and  57.   All  of  these 
components  are  at  odd  sequencies,  so  the  power  will  be  located 
in  the  odd  number  sequencies  of  the  PSD,  and  the  odd  numbered 
components  will  combine  into  one  component  of  the  Group 
Spectrum.   Again,  see  Equations  (3.34). 

The  FWT,  PSD,  and  Group  Spectrum  of  eight 
cycles  of  a  sine  wave  in  the  Walsh  interval  is  analyzed  in 
Figures  3.14  to  3.16.  The  FWT  is  considerably  simpler,  with 
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Figure  3.9.   PSD  o-f  1  Cycle  o-f  a  Sine  Function. 
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Group   Spectrum   o-f     1    Cycle   o-f    a    Sine   Function 
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Figure    3.11.       FWT   o-f    3   Cycles   o-f    a   Sine    Function 
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Figure   3.12.       PSD   o-f    3   Cycles   o-f    a   Sine   Function 
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Figure   3.13.       Group   Spectrum   o-f    3   Cycles   o-f    a   Sine   Function 
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Figure    3.14.       FWT   o-f    8   Cycles   o-f    a   Sine    Function. 
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Figure   3.15.       PSD    o-f    8   Cycles    o-f    a   Sine   Function. 


97 


I— 
CD 
LU 
Q, 
CO 

0_ 
ZD 

o 

Ci 
CD 


CO 


a 

CO 

Hi 

_J 
o 

>- 

CJ 
GO 


1  CD 


in 


00 


CM 


Ui 

CD 


a 

CD 


Si 


I L 


J i i L 


'-•ocDc^cDin^cncM 


S3 


CSZ8I  "=D    OVW   ISA 


Figure  3.16.   Group  Spectrum  o-f  8  Cycles  o-f  a  Sine  Function 
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only  4  components  at  PSD  sequencies  of  8  and  24.   The  PSD 
shows  the  peak  power  at  a  sequency  o-f  8  and  a  notable  harmonic 
at  24.  The  Group  Spectrum  -follows  -from  Equations  3.34. 

A  few  conclusions  about  the  transform  of 
a  sinusoid  could  be  drawn.   The  peak  coefficient  of  the  FWT 
falls  at  a  Walsh  index  that  translates  to  a  sequency  that  is 
equal  to  the  sinusoidal  frequency.   The  PSD  then  has  a 
component  that  is  maximum  at  this  sequency.  In  addition,  this 
maximum  power  component  is  usually  accompanied  by  notable 
components  at  harmonics  of  this  sequency. 

The  Group  Spectrum,  although  highly 
compressed  and  uncomplicated,  requires  a  knowledge  of  the  PSD 
for  it  to  be  useful .   Note  that  the  Group  Spectrum  of  the  1 
cycle  sine  wave  is  equal  to  the  the  Group  Spectrum  of  the  3 
cycle  wave.    An  educated  guess  says  that  the  Group  Spectrum 
for  all  odd  numbers  of  cycles  would  be  the  same,  that  of  Gj , 
because  only  odd  numbered  components  <Eq.  3.34)  of  sequency 
are  present.   Each  component  of  the  Group  Spectrum  could 
represent  a  series  of  sinusoidal  frequencies,  with  no  way  to 
distinguish  the  actual  frequency  of  the  sinusoid  without  prior 
knowledge  of  the  PSD. 

b.   The  FWT  of  a  Rectangular  Function. 

Consider  the  FWT,  PSD,  and  Group  Spectrum  of  a 
4  hertz  square  wave.   Knowing  that  the  Walsh  function  lends 
itself  readily  to  the  reproduction  of  such  a  wave,  simple 
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results  -for  graphs  o-f  the  trans-form  and  related  computations 
would  be  expected. 

Indeed  this  is  the  case.   Figures  3.17  and 
3.18  show  the  simplicity  o-f  representation  -for  rectangular 
-functions.   The  FWT  component  at  n  =  7,  sequency  =  4,  is  the 
only  coe-f-f icient  in  Figure  3.17.   (Note  that  a  4  hertz  square 
wave  is  WAL<7,t)).   All  o-f  the  power  lies  at  a  sequency  o-f  4 
(Figure  3.18),  and  in  a  Group  Number  3  (Eqn.  3.34).   With  a 
proper  number  o-f  samples,  any  square  wave  of  a  particular 
■frequency  would  be  represented  with  similar  characteristics. 

A  rectangular  -function,  page  73,  composed  o-f 
square  pulses  was  expanded  in  a  Walsh  Series  in  this  chapter. 
Consider  the  FWT,  PSD,  and  Group  Spectrum  o-f  the  -function. 

Figures  3.19,  3.28  and  3.21  show  the  results 
o-f  the  computation.  These  graphs  are  not  normalized  in 
order  to  show  the  similarities  with  the  previously  computed 
example  on  page  73.    The  FWT  shows  2  coe-f -f icien ts  o-f  8.5 
magnitude,  at  n  =  8  and  n  =  2,   matching  the  results 
determined  be-fore.   Since  the  FWT  outputs  coe-f  -f  icien  ts  based 
on  positive  phasing,  this  result  is  correct. 

The  power  contained  in  the  function  is  in  the 
DC  term  and  the  first  sequency  component  as  shown  in  Figure 
3.28.   These  sequency  components  result  in  2  components  in  the 
Group  Spectrum.   See  Equations  3.34. 
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Figure  3.17.   FWT  o-f  a  4  Cycle  Square  Nave 
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Figure  3.18.   PSD  o-f  a  4  Cycle  Square  Wave 
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Figure    3.28.       PSD    o-f    2    Rectangular    Pulses 
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Figure  3.21.   Group  Spectrum  o-f  2  Rectangular  Pulses 
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c.   The  FWT  and  a  Input  Time  Shi-ft. 

Section  4-c  stated  that  the  Discrete  Walsh 
Trans-form  is  not  invariant  to  a  time  shi-ft  or  phase  shift  o-f 
the  input.   Let  sin<t)  be  circularly  shi-fted  K/4  to 
sin<t-*/4).   Now  examine  Figures  3.22,  3.23,  and  3.24  and 
compare  with  Figures  3.8,  3.9,  and  3.18. 

Note  the  shi-fted  FWT  coe-f -f i cien ts  pattern  is 
more  complex.   A  sizable  component  now  lies  at  n  =  2.   This  is 
a  CAL<l,t)  component  with  a  sequency  o-f  1.   The  unshi-fted  FWT 
graph  contained  only  a  negligible  CAL<  1 ,  t)  component. 
Double  components  at  n  =  5  and  6,  13  and  14,  and  29  and  39 
represent  sequencies  o-f  3,  7,  and  15,  respectively,  the  same 
components  that  were  contained  in  the  unshi-fted  sinusoid,  but 
the  magnitudes  o-f  the  components  at  the  same  index  Cn3  are  not 
equal . 

Although  the  FWT  magnitudes  are  di-f-ferent  at 
each  index  Cn3  o-f  the  shifted  and  unshifted  sinusoids, 
remember  that  adjacent  components  combine  to  form  a  PSD 
coefficient  of  a  particular  sequency.   The  components  vary  in 
a  reciprocal  manner,  so  when  they  are  squared  and  summed,  the 
time  variant  effect  is  indeed  lessened  [23:  p. 89] .   The  result 
is  a  PSD  that  has  sequency  components  that  are  \/ery   close  to 
being  time  invariant. 

One  expectation  that  was  revealed  to  be  true 
was  the  time  invariance  of  the  Group  Spectrum.   Beauchamp  C23: 
p.  186]  states  this  spectrum  is  time  invariant  to  time  shifts 
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Figure    3.22.       FWT   o-f    a   Shi-fted    Sinusoid 
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Figure  3.23.   PSD  o-f  a  Shi-fted  Sinusoid 
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Figure  3.24.   Group  Spectrum  o-f  a  Shifted  Sinusoid 
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of  the  input,  and  examination  of  Figure  3.23  and  its  magnitude 
gives  an  a-f-f irmative  response  to  his  assertion. 

The  rectangular  pulse  -function  will  positively 
demonstrate  the  effects  o-f  time  shift.   The  FWT,  PSD,  and 
Group  Spectrum  of  Figure  <3.25a)  is  shown  in  Figures  3.19, 
3.29,  and  3.21.   The  FWT  and  associated  spectrums  of  Figure 
(3.25b)  is  shown  in  3.26,  3.27,  and  3.28. 

Quite  a  difference  in  the  FWT  of  the  shifted 
function  can  be  noticed.   The  coefficient  pattern  is  much  more 
"spread  out",  and  the  magnitudes  of  the  prominent  coefficients 
at  n  =  8  and  1  are  less  in  the  shifted  rectangular  function. 

The  PSD's  are  different  also.  With  the  power 
being  spread  among  more  sequencies  in  the  shifted  rectangular 
function,  the  magnitudes  of  the  power  coefficients  at  8  and  1 
are  reduced. 

The  Group  Spectrums  are  not  alike  in  this 
case.   The  two  input  sequences  to  the  FWT  program  are  made  up 
of  samples  of  the  original  function  and  the  shifted  function. 
These  sequences  are  not  alike  and  generate  different  Group 
Spectrums. 

Has  Beauchamp's  assertion  C23:  p. 106]  been 
violated?  The  answer  is  no.   In  the  sinusoidal  shift,  it  was 
assumed  that  the  function  was  periodic  outside  the  interval 
8<t<l  over  which  the  Walsh  functions  are  defined  and  over 
which  the  function  was  taken.   When  the  sinusoid  was  shifted, 
the  part  shifted  out  of  the  interval  "wrapped  around"  into  the 
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Figure  3.25.   Two  Rectangular  Pulses,  Shifted  and  Unshi-fted 
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Figure  3.26.   FWT  o-f  Shifted  Rectangular  Pulses 
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Figure  3.27.   PSD  o-f  Shifted  Rectangular  Pulses 


113 


CD 


\m 


en 
u. 

LL. 
LU 

a 
u 

u 
m 

a. 


CD 


CD 


I— 
U 
UJ 


en 


♦  CO 

I       a: 

I       ca 

i       =3 

I  CM    Z 

a. 

ZD 

a 

a: 

CD 


4  S3 


*-* 

CD 

CD 

t** 

CO 

in 

• 

S3 

IS 

s 

S3 

S3 

S3 


CO 

S3 


CM 
S3 


S3 


S3 


ovw 


Figure  3.28.   Group  Spectrum  o-f  Shi-fted  Rectangular  Pulses 
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■first  part  o-f  the  interval  (circular  shift)  ,  resulting  in  the 
shi-fted  -function  having  the  same  power  as  the  unshi-fted 
•function  . 

However,  the  rectangular  -function  wasn't 
de-fined  as  a  periodic  -function,  merely  two  pulses  in  an 
interval.   When  the  rectangular  function  was  shi-fted,  the  part 
that  shi-fted  outside  the  interval  wasn't  wrapped  around  to  the 
■front  (not  circularly  shi-fted)  and  was  lost,  decreasing  the 
total  power  found  in  the  shifted  function.  The  PSD  or  the 
Group  Spectrum  could  not  be  the  same  as  the  unshifted 
function,  since  the  amount  of  total  power  is  different  in  each 
of  the  functions. 

A  conclusion  is  made  that  the  Group  Spectrum 
is  time  invariant  to  phase  shifts  in  the  input  for  periodic 
circular  shifted  functions  only. 
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IV.   THE  HALSH  TRANSFORM  AND  DE INTERLEAVING 

The  -first  two  chapters  presented  some  background  for  the 
appreciation  of  the  role  of  ESM  and  the  basic  components  o-f  a 
typical  ESM  system.  Understanding  the  tools  and  purpose  of  the 
trade,  though,  only  gives  a  small  insight  into  the  nature  of 
probl em. 

There's  a  lot  of  electromagnetic  radiation  out  there!   The 
combat  environment  presents  an  overwhelming  amount  of  energy 
that  has  to  be  collected,  analyzed,  and  sorted.   The  data 
stream  pours  into  the  system  from  all  angles,  with  numerous 
frequencies  and  amplitudes.   Corruptions,  disturbances, 
reflections,  and  missing  pulses  make  the  dein ter 1 eaving 
process  a  most  difficult  problem.  Today's  techniques, 
described  in  Chapter  Two,  do  a  good  job,  but  improvements  are 
needed  and  sought  after. 

Unique  in  many  respects,  the  Nalsh  Functions,  which  are 
described  in  Chapter  Three,  are  examined  for  their  usefulness 
in  deinterl eaving  pulse  trains  into  separate  chains  of  pulses 
with  different  pulse  repetition  intervals.   More  specifically, 
this  chapter  attempts  to  determine  if  there  is  a  PRI/PRF 
recognition  feature  in  the  Mai sh  Transform  of  the  data 
presented  by  an  ESM  receiver. 
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A.   REPRESENTATION  OF  ESM  RECEIVER  DATA 

The  operating  element  o-f  the  typical  deinterl eavi ng 
process  is  a  computer  data  word  (or  words)  o-f  some  -format 
appropriate  to  the  system. 

The  data  word  contains  values  o-f  the  initial  parameters 
measured  by  the  system,  and  basically  these  values  are 
examined  for  matches  in  previously  analyzed  pulses.  So-ftware 
processing  systems  make  these  comparisons  and  places  words  o-f 
similar  parameters  into  cells  or  histograms,  patiently  waiting 
until  enough  pulse  information  can  be  -found  to  make  an 
identification.    Other  algorithms  process  the  data  word  and 
make  calculations  of  PRI's,  adding  a  new  sorting  parameter 
that  is  used  in  the  dein ter 1 eaving  process. 

The  investigation  of  deinterl eaving  with  the  FWT  will  need 
a  different  symbol ogy  for  the  receiver  data,  since  the  FWT 
operates  on  a  series  of  N  numbers,  not  a  data  word,  to  obtain 
the  transform. 

Imagine  an  absolute  time  line  that  consists  of  some  point 
in  the  past,  and  ends  at  the  present.   Consider  a  received 
pulse  as  being  represented  by  a  vertical  line  stationed  at 
the  time  of  arrival  (TOA)  of  that  pulse.   The  height  of  the 
line  could  be  made  proportional  to  the  amplitude  of  the 
received  pulse.   The  point  on  the  time  line  where  the  pulse 
•stands'  is  the  time  of  arrival  of  the  pulse. 

Figure  4.1  looks  like  a  messy  situation,  but  there  are 
really  only  four  pulse  trains  of  different  PRI's  on  the  line. 
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Figure  4.1.   TOA  representation  o-f  received  pulses. 

PRI's  could  be  calculated  by  subtracting  successive  TOA's 
o-f  pulses  in  the  same  pulse  train. 

PRI  !  =  T0A12  ~  TOAn 

PRI2  =  T0A22  "  T0A21 
PRI 3  =  T0A32  -  T0A31 
PRI 4  -  T0A42  -  T0A41 

This  procedure  is  used  when  other  parameters  are  present 
to  identi-fy  the  pulses. 

Keep  in  mind  that  there  is  a  steady  stream  o-f  these  pulses 
into  the  ESM  system.   In  order  to  make  things  manageable,  the 
number  o-f  pulses  being  examined  has  to  be  reduced. 

Leml ey  C28  :  p. 12]  calls  this  reduction  mechanism  the 
"presorting  aperture."   It  is  o-f  variable  width  to  accomodate 
high  and  low  density  environments  and  is  similar  to  the  time 
slice  concept  discussed  under  deinter 1 eaving  in  Chapter  Two. 
It  reduces  the  number  o-f  pulses  being  processed  at  one  time. 

Another  way  of  thinking  o-f  this  reduction  concept  is  to 
visualize  a  "window"  on  the  stream  o-f  pulse  lines.   Only  a 
certain  number  o-f  pulse  lines  can  be  present  in  the  window  at 
anytime,  and  deinter 1 eaving  operations  are  carried  out  on 
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these  pulse  lnes.   The  window  can  be  divided  into  small  time 
intervals,  which  could  be  analogous  to  the  TOA  resolution 
specification  o-f  a  receiver  (usually  198-288  nanoseconds).  I-f 
a  pulse  line  is  present  in  the  time  interval,  a  number 
proportional  to  the  amplitude  is  generated  for  the  pulse  line 
representation.   If  not,  a  zero  is  used. 

The  author  prefers  to  visualize  that  the  receipt  of  a 
pulse  triggers  the  collection  of  the  next  N  time  intervals  -for 
analysis  and  deinter 1 eaving.   Each  interval  could  contain  a 
pulse  or  not. 

B.   CONSTANT  AMPLITUDE  TOA  STRINGS  AND  THE  FWT 

This  chapter  will  work  mostly  with  N  =  64  intervals  and 
pulse  line  representations.  N  must  be  a  power  of  two  to  be 
used  in  the  FWT. 

The  pulse  line  representation  can  be  symbolized  for 
computer  work  with  the  BASIC  language  DATA  statement.   Using 
the  statement  on  the  64  pulse  lines  in  the  window  gives  a 
representation  that  accounts  for  the  two  variables  that  are 
present  in  the  pulse  line  stream,  the  amplitude  and  TOA  of  the 
pul ses. 

Such  a  DATA  line  for  Figure  4.1  would  be 

vt-i  In    1  ,  .  O  |  .  'd  ,  .  0)t)  |  •  D  ,  ■  /  o  ,  1,0,  •  o  ,  •  O  ,  .0,0,  >  O  ,  1,  .5,0,  •  O  ,  .  f  5 

<4.  1) 
Each  number  represents  the  amplitude  of  a  recognized  pulse 
that  occupies  a  TOA  time  interval . 
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Think  o-f  it  in  one  o-f  two  ways: 

1.  The  receiver  outputs  a  measured  value  o-f  a  pulse 
amplitude  each  TOA  resolution  time  interval.   I-f  no  pulse  was 
present,  it  outputs  a  zero.   Or, 

2.  The  receiver  outputs  a  pulse  amplitude  only  upon 
detection  and  analysis  o-f  a  pulse.   Time  intervals  o-f  TOA 
resolution  width  in  between  pulse  detections  can  be 
represented  with  zeroes. 

Either  way,  the  DATA  line  looks  like  (4.1). 

Assume  -for  simulation's  sake  that  the  period  between  DATA 
numbers  is  a  known  value  o-f  time,  and  the  DATA  line  represents 
the  output  o-f  an  ESM  receiver  to  the  preprocessor.   The 
numbers  indicate  amplitude  o-f  the  received  pulses,  and  their 
distance  in  intervals  between  similar  pulses  is  a  measure  o-f 
their  PRI's. 

Things  can  be  simplified  a  bit  -further  by  considering  all 
pulse  lines  o-f  the  TOA  string  to  be  o-f  the  same  amplitude. 
Only  one  variable  would  then  be  present  in  the  string,  the 
TOA. 

A  DATA  line  representation  o-f  a  string  o-f  this  type  is 

DATA   1,1,1,1,9,1,1,1,8,1,1,1,8,1,1,1,8,1,1,  (4.2) 

In  short,  this  TOA  string  and  its  DATA  representation  is 
merely  an  indication  o-f  the  receipt  o-f  a  pulse  or  not  during 
the  resolution  time  interval. 

A  number  o-f  Fast  Walsh  Trans-forms  were  run  on  simulated 
DATA  TOA  strings,  and  plots  o-f  the  coe-f -f icients  were  made.  The 
Power  Spectral  Density  and  Group  Spectrum  coe-f-f icients  were 
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also  calculated  and  plotted.  All  o-f  the  Group  Spectrum 
plots  are  not  included,  as  the  emphasis  o-f  this  examination 
was  placed  on  the  FWT  and  PSD.   Some  o-f  them  are  discussed  in 
Section  C  o-f  this  chapter. 

The  plots  are  grouped  by  type  <FWT  or  PSD)  in  Appendix 
D.  Each  plot  is  normalized  by  dividing  all  coe-f-f icien ts  by 
the  maximum  coe-f -f  icien  t .  Actual  values  o-f  the  coe-f-f  icien  ts 
can  be  determined  with  the  scale  on  the  vertical  axis.  The 
maximum  component's  magnitude  is  shown  in  parenthesis. 

These  plots  show  the  index  number  and  sequency  o-f  the  FWT 
and  PSD  components,  respectively,  but  only  their  relative 
magnitude  with  respect  to  the  maximum  component.  Keep  that  in 
mind  when  examining  the  plots.   Where  important,  the  FWT  and 
PSD  are  replotted  with  the  actual  magnitudes  on  the  vertical 
axis. 

These  DATA  representations  simulate  single  PRI  pulse 
trains  received  by  an  ESM  system.   Note  the  number  o-f 
intervals  between  successive  pulses  is  indicated  in  the  title 
o-f  the  plot.   To  calculate  a  simulated  PRI  o-f  the  train, 
simply  multiply  the  number  o-f  intervals  between  sucessive 
pulses  by  the  time  interval  resolution,  a  known  value  and  a 
■function  o-f  the  receiver. 

The  plots  are  examined  by  type  in  the  -following  sections, 
and  conclusions  are  reached  and  reported.   These  single  PRI 
representations  are  the  building  blocks  for  the  interleaved 
pulses  that  are  investigated  later  in  the  chapter. 
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1.   Fast  Walsh  Trans-forms  o-f  the  TOA  Strings 

An  examination  o-f  the  FWT's  separated  them  into  two 
groups  of  similar  characteristics. 

The  -first  group  is  composed  o-f  FWT's  with  intervals 
between  pulses  that  are  even  numbers.   Examine  each  plot 
in  Appendix  D  (page  176  to  197)  and  note  that  each  even 
interval  FWT  is  symmetrically  even  about  a  point  midway 
between  n  =  31  and  32.   Also,  i-f  the  even  number  is  a  power  o-f 
two  (p.o.t.),  the  FWT  is  considerably  less  complex.   Each 
successive  p.o.t.  interval  is  composed  of  the  addition  of  a 
pair  or  pair  of  coefficients  to  the  previous  p.o.t.  FWT,  and 
its  magnitude  is  one  half  the  previous  p.o.t.  FWT. 

For  example,  one  pair  of  coefficients  is  added  to  TOA 
TAG  2  at  n  =  31  and  32  (sequency  of  16)  to  form  TOA  TAG  4. 
Two  pairs  of  coefficients,  (n  =  15  and  16,  n  =  47  and  48)  are 
added  to  TOA  TAG  4  to  make  TOA  TAG  8,  4  pairs  to  make  TOA  TAG 
16,  etc.   At  N  —  64,  enough  pairs  of  coefficients  have  been 
added  to  have  one  coefficient  for  each  En]. 

The  odd  numbered  intervals  do  not  have  symmetrical 
characteristics.   In  addition,  the  highest  relative  magnitude 
Walsh  index  component  or  components  seem  to  bear  no  relation 
to  the  interval  number  and  thus  the  PRI  of  a  pulse  train.   <An 
unfortunate  conclusion  for  PRI  recognition.) 

Close  scrutiny  of  both  appearance  groupings  was  made 
in  an  effort  to  notice  a  feature  that  is  distinctive  to  the 
particular  similated  PRI.   There  doesn't  seem  to  be  any 
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present,  either  in  magnitude  or  coe-f-f icient  distribution,  that 
is  obvious  without  a  detailed  examination. 

0-f  initial  interest  to  the  author  was  the  magnitude  of 
the  WAL<8,t)  component.   This  magnitude  seemed  to  be  related 
to  the  number  o-f  resolution  intervals  between  the  pulses.   It 
holds  -for  most  but  not  all  cases.   Consider  Table  4.1  which 
lists  the  TAG  number  and  the  WAL<8,t)  component. 

Table  4.1.   Amplitude  o-f  WAL<8,t>  and  TAG  No.  Inverse 
TOA  TAG  no.         HAL<8,t)  amplitude         TAG  no.  inverse 


1 

1.8 

1.8 

2 

8.5 

8.5 

3 

8.34375 

8.3333 

4 

8.25 

8.25 

5 

8.283125 

8.28 

6 

8. 171875 

8.  1667 

7 

8. 15625 

8. 1429 

3 

8.  125 

8.  125 

9 

8.  125 

8.1111 

18 

8.  189375 

8.  1 

11 

8.89375 

8.8989 

12 

8.89375 

8.8833 

13 

8.878125 

8.876923 

14 

8.878125 

8.871429 

15 

8.878125 

8.8667 

16 

8.8625 

8.8625 

19 

8.8625 

8.8526 

22 

8.846875 

8.845455 

25 

8.846875 

8.84 

32 

8.83125 

8.83125 

48 

8.83125 

8.825 

51 

8.83125 

8.8  19688 

64 

8.8  15625 

8.8  15625 

For  instance,  TOA  TAG  2  has  two  intervals  between 
pulses,  and  the  magnitudeo-f  WAL(8,t)  is  8.5,  or  the  inverse  o-f 
2.   TOA  TAG  3  has  a  magnitude  o-f  .34375,  which  is  close  to  the 
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inverse  o-f  3.   TOA  TAG  4  has  a  WAL<8,t)  amplitude  of  .25, 
which  is  one  divided  by  -four,  the  number  o-f  intervals  between 
pulses  in  the  representation. 

I-f  there  was  a  relationship  between  the  number  o-f 
intervals  between  the  pulses  in  the  DATA  line  and  the 
amplitude,  then  perhaps  this  would  carry  over  into  the 
interleaved  pulse  trains  composed  of  these  single  PRI 
representations.   Exceptions  are  evident,  though,  such  as  TOA 
TAG  9,  which  has  the  same  magnitude  as  TOA  TAG  8,  although  the 
coefficient  distribution  is  not  the  same.  This  implies  that 
any  recognition  feature  will  probably  depend  upon  two 
variables,  perhaps  the  amplitude  and  the  sequency  of  the 
coefficients. 

2.   Power  Spectral  Densities  of  the  TOA  Strings 

An  examination  of  the  PSD's  in  Appendix  D,  page 
198  through  page  219,  will  yield  similar  conclusions. 

PSD's  of  even  interval ed  TOA  strings  are  symmetrically 
even  about  a  sequency  of  16.       PSD's  of  even  numbers  that  are 
also  a  p.o.t.  are  considerable  less  complex,  and  are  one 
fourth  the  magnitude  of  the  previous  p.o.t.  PSD  (except  for 
TOA  TAG  2  to  TOA  TAG  4) . 

The  odd  numbered  PSD's  have  no  symmetry,  and  sequency 
components  that  stand  out  in  relative  magnitude  do  not  bear  a 
relationship  to  the  simulated  PRI  of  the  string.   (Remember 
the  number  of  intervals  between  the  TAG's  is  a  measure  of  the 
PRI  of  the  simulated  pulse  train,  when  the  time  interval 
between  numbers  is  known.) 
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3.      General    Comments 

The   entire   purpose   of    this   examination    has   been    to 
recognize   similarities   or    relationships   between    succcessive 
TOA  TAG   representations    in    order    to   recognize    the   same 
•features    (if    they   exist)     in    an    interleaved   pulse    train    DATA 
line   representation. 

At    this   point,    the   results   are  very   general    and 
an    optimistic   mood   regarding   any   usable   -features  might    belong 
only    to    the   optimist.       It    seems    that    the   even    numbered    graphs 
<FWT   and   PSD   )    are    the   only   graphs    that    show   any   usable 
■features    that    could   be   used   -for    recogniton.      One   shouldn't 
give   up   all    hope,    though. 

The   author    believes    it    to   be   possible    to   develop   a 
computer    algorithm    that   would   recognize   and   determine    the  TOA 
TAG   number    (and    thus   a   PRI    o-f    a   simulated   pulse    train)    for    a 
given    single   pulse    train    DATA    line   representation.      Enough 
features   exist    in    the   FWT   and   PSD    that    could   be   used    to 
process   a   particular    single   PRI    representation   with    a   Known 
resolution    time    interval    and   output    the   correct    PRI.    The 
program  would   have    to   determine    the   PRI    by  matching 
characteristics    (sequency   components   and  magnitudes)    of    the 
FWT   and/or    PSD    to   a   data   bank    composed   of    characteristics   of 
the   FWT's   and   PSD's   of    single   pulse    train    representations. 
For    example,    it    could   first    take    into   account    the   symmetry   of 
the   FWT/PSD.       If    it    is   symmetrically   even,    then    it    isn't    an 
odd    interval    TOA  TAG. 
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The  algorithm  would  be  reasonably  complex  -for  just 
this  simple  TOA  TAG  representation  of  a  single  pulse  train  of 
a  particular  PRI .   But  it  could  be  done. 

Since  it  addresses  only  single  PRI  TAG 
representations,  the  algorithm  would  only  be  effective  if  you 
could  convince  an  opponent  to  use  only  one  radar  at  a  time 
against  the  ESM  platform.   Obviously,  this  isn't  practical. 

More  realistically,  an  examination  o-f  interleaved 
pulse  trains  and  a  determination  o-f  a  PRI  recognition  -feature, 
if  one  exists,   would  be  in  order.   This  will  be  done  in 
Section  D.  after  discussion  of  the  effects  of  a  time  shift  in 
the  TOA  TAG  string. 

Some  limitations  of  the  DATA  line  symbol ogy  should  be 
stated.   Using  N  =  64  and  a  resolution  between  numbers  of  208 
nanoseconds  means  the  window  length  is  less  than  2 
microseconds  long.   This  is  not  a  realistic  figure,  since  a 
high  density  of  19*  pulses  per  second  gives  an  average  of  1 
pulse  every  microsecond.   The  use  of  N  =  1024  for  an 
interleaved  DATA  line  representation  might  be  far  better  in 
the  numbers  of  coefficients  <and  more  realistic  in  length  of 
the  window),  but  unless  an  identifiable  PRI  feature  is 
present,  the  extra  data  is  exactly  that,  extra  data. 

An  additional  limitation  of  this  representation 
concerns  coincident  pulses.   Pulses  arriving  at  the  receiver 
at  the  same  time  (falling  in  the  same  resolution  time 
interval)  are  not  individually  represented  in  the  DATA  line 
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symbol ogy.   For  this  reason,  one  cannot  represent  a  TOA  TAG  2 
and  any  other  even  numbered  TAG.  This  limitation  is  a 
problem  that  will  affect  the  FWT  and  PSD  o-F  the  interleaved 
pulse  train  representation. 

C.   TIME  SHIFTED  TOA  STRINGS 

Chapter  Three  has  stated  that  the  FWT  and  PSD  are  not  time 
invariant  to  circular  shifts  o-f  the  input  number  series. 
Figures  4.2  through  4.5  are  FWT's  o-f  shifted  TOA  TAGS  . 
Compare  them  with  Figures  D.4,  D.ll,  D.13,  and  D.16. 

The  magnitudes  o-f  the  FWT  coe-f-f icients  at  each  Cn]  are  not 
the  same  in  the  shifted  TAGS.   Compare  TOA  TAG  5  and  TOA  TAG  5 
SHF  1.   However,  notice  that  shifted  TOA  TAG  12  <TOA  TAG  12 
SHF  1)  has  the  same  magnitude  at  each  WALSH  index  Cn]  as  does 
the  unshifted  TOA  TAG,  but  the  components  from  n  -  32  to  n  = 
63  are  opposite  in  sign.  The  even  symmetry  of  the  even 
numbered  FWT  about  the  midpoint  has  been  changed  to  odd 
symmetry  about  midpoint. 

This  is  a  \fery    neat  and  clean  distinction,  and  further 
emphasizes  the  distinction  between  even  and  odd  numbered  TOA 
TAGS.   It  is  further  evidence  that  at  least  some  general 
features  exist  which  could  be  made  to  work  in  a  computer 
algorithm  to  identify  individual  PRI's. 

FWT's  of  odd  numbered  intervals  have  no  general 
distinctive  change  between  shifted  plots  and  standard  plots. 
In  a  few  shifted  odd  numbered  plots,  there  are  some  particular 
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Figure  4.2.   FWT  :  TOA  TAG  5  SHIFT  1 
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Figure  4.3.   FWT  :  TOA  TAG  12  SHIFT  1 


12? 


Li- 

t—i 

en 

z: 

LL 

< 

T* 

C£ 

CO 

t— 

^J- 

X 

r-t 

en 

_j 

CD 

<: 

< 

3c 

h- 

H- 

< 

CO 

a 

< 

h- 

Ll. 

in 


ca 


in 

« 
1 


-CD 


CD 

in 


GO 


CM 
CO 


CD 


3 


CSI8Z0'-I)    SJVW  13d 


Figure  4.4.   FWT  :  TOA  TAG  14  SHIFT  1 
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Figure    4.5.       FWT    :    TOA   TAG    19    SHIFT    1 


131 


magnitudes  of  index  Cn]  that  remain  the  same  or  increase 
proportionately,  but  others  do  not.   If  any  real  relationship 
between  the  shifted  and  unshifted  odd  numbered  FWT  plots 
exist,  then  it  will  have  to  be  revealed  with  a  much  more 
detailed  study. 

Quite  simply,  odd  numbered  shi-fted  FWT  plots  show  no  major 
distinctive  features  that  are  readily  apparent.   Even  numbered 
shifted  FWT  plots  follow  the  sign  change  rule  for  components 
at  n  =  32  to  63,  but  the  magnitudes  of  the  components 
of  both  plots  are  the  same. 

Although  the  odd  numbered  PSD  plots  change  in  the 
amplitude  of  the  sequency  components  from  unshifted  to 
shifted,  the  even  numbered  PSD's  are  the  same  in  both  shifted 
and  unshifted  plots.   Indeed,  they  are  the  same  in  both 
magnitude  and  location  of  sequency  components.   Compare 
Figures  4.6  to  4.9  on  the  following  pages  with  Figures  D.26, 
D.33,  D.35,  and  D.38  in  Appendix  D. 

Thus,  the  even  numbered  PSD's  of  this  DATA  line 
representation  join  the  time  invariant  Group  Spectrums 
in  being  immune  to  the  effects  of  a  circular  time  shift. 
Compare  the  shifted  Group  Spectrums  Figures  4.18  to  4.13,  with 
the  unshifted  Group  Spectrums  in  Appendix  E,  Figures  E. 1 
through  E.4. 

Again  the  point  is  made  that  these  general  features  could 
be  recognized  in  a  PRI  recognition  algorithm.   At  the  \^ery 
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Figure  4.6.   PSD  :  TOA  TAG  5  SHIFT  1 
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Figure  4.7.   PSD  :  TOA  TAG  12  SHIFT  1 
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Figure  4.8.   PSD  :  TOA  TAG  14  SHIFT  1 
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Figure  4.9.   PSD  :  TOA  TAG  1?  SHIFT  1 
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Figure  4.18.   GROUP  SPECTRUM  :  TOA  TAG  5  SHIFT  1 
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Figure  4.11.   GROUP  SPECTRUM  :  TOA  TAG  12  SHIFT  1 
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Figure  4.12.   GROUP  SPECTRUM  :  TOA  TAG  14  SHIFT  1 
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Figure    4.13.       GROUP    SPECTRUM    :    TOA   TAG    1?    SHIFT    1 
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least,  it  would  eliminate  the  odd  or  even  choices  of  time 
intervals  between  numbers  of  the  single  PRI  representation. 

D.   INTERLEAVED  TOA  STRINGS  AND  THE  FWT 

The  receipt  o-f  radar  pulses  from  several  di-f-ferent 
emitters  results  in  individual  pulse  trains  that  are 
interleaved  into  a  data  stream  -from  the  receiver.   Separation 
or  dein ter  1  eaving  of  these  pulse  trains  -from  the  data  stream 
must  be  done  to  identi-fy  the  individual  emitters.    This 
section  examines  the  plots  o-f  the  FWT  and  PSD  coe-f-f icients  o-f 
a  DATA  line  representation  o-f  interleaved  pulse  trains 
composed  o-f  the  TOA  TAG  representation  in  Section  B.  o-f  this 
chapter . 

It  is  desired  to  answer  the  -following  questions: 

1.  Are  the  components  o-f  the  FWT  and  PSD  o-f  an 
interleaved  representation  identifiable  as  belonging  to  a 
particular  PRI? 

2.  Is  there  a  relationship  between  the  magnitude  or 
sequency  o-f  the  components  o-f  an  interleaved  plot  with  that 
o-f  plots  o-f  single  PRI  trains? 

3.  Are  interleaved  plots  made  by  simply  adding  single  PRI 
plots? 

Begin  by  considering  interleaved  TOA  TAGS  o-f  Table  4.2. 

The  DATA  line  representation  -for  TOA  TAG  5  wi  th  N  =  64  is 

DATA  1,8,9,8,9,1,8,8,8,8,1,8,8,8,9,1,8,9,8,8,1,3,9,9,8,1,8, 
9,8,8,  1,9,9,9,9,  1,9,9,9,9,  1,9,9,9,9,  1,9,9,9,9,  1,8,9,9, 
9,1,9,9,9,9,1,9,9,9  (4.3) 

and  for  TOA  TAG  8, 

DATA  1,9,9,9,9,9,9,9,1,9,9,9,9,9,9,9,1,9,9,9,9,9,9,9,1,9,9, 
9,9,9,9,9,1,9,9,3,8,8,8,8,1,8,8,8,8,8,8,8,1,8,8,8,8,8, 
8,3,1,8,3,8,3,8,3,8  (4.4) 
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Now,  interleave  the  two, 

DATA  1,8,8,0,8,1,8,8,1,8,1,8,9,8,9,1,1,8,8,8,1,8,0,9,1,1,8, 
8,8,8,1,8,1,8,8,1,8,8,8,8,1,8,8,8,8,1,8,8,1,8,1,8,9,0, 
8,1,1,8,8,8,1,8,8,8  <  4 . 5) 

Notice  that  the  -first  pulse  of  each  TAG  is  coincident  and 
both  are  represented  by  the  -first  pulse  o-f  the  interleaved 
DATA  1 ine. 

The  FWT  o-f  TAGS  5  and  3  are  shown  in  Figures  4.14  and 
D.7.   TOA  TAG  5  components  are  emphasized  by  setting  the 
WAL(8,t>  components  equal  to  8  be-fore  plotting.   The  plotting 
routine  then  normalizes  to  a  smaller  magnitude  component  and 
increases  the  relative  magnitude  o-f  the  other  components. 

Examine  the  interleaved  FWT,  Figure  4.15.   Note  sizable 
components  at  9,  22,  25,  31,  32,  37,  45,  48,  58,  68,  and  63. 
Now  examine  TOA  TAG  5  with  WAL(8,t)  equal  to  8.   Note 
emphasized  components  at  25,  37,  44,  45,  58,  51,  and  52. 

There's  some  commonality  between  the  plots.   Components  at 
25,  37,  45,  58,  and  52  are  obviously  associated  with  TOA  TAG 
5.   So  where  is  TOA  TAG  8?  A  closer  look  reveals  its 
in-fluence  on  the  interleaved  FWT.   Note  identical  relative 
magnitude  pairs  at  n  =  15  and  16  and  n  =  31  and  32.   These  are 
generated  by  the  coefficient  pairs  of  these  [n]  from  TOA  TAG 
8.   Note  also  that  n  =  48  of  the  interleaved  FWT  is  higher  in 
relative  magnitude  because  there  is  a  small  n  =  48  component 
of  TOA  TAG  5. 

How  would  one  use  this  interleaved  FWT  to  deinter leave 
pulse  trains?  A  direct  computation  method  using  a  correlation 
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Figure.    4.14.       FWT    :    TOA  TAG    5 
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Figure.  4.15.   FWT  :  TOA  TAG  5  AND  8  INTERLEAVED 
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algorithm  is  considered  in  the  next  section.   The  author 
suggests  a  "matching-  algorithm  could  be  written  that  would 
scan  the  interleaved  FWT,  and  note  relatively  high  magnitudes 
at  certain  Cn]  above  some  threshold.   These  Cn]  could  be 
matched  to  a  data  bank  of  important  Cn]  components  of  single 
PRI  FWT  plots,  along  with  any  special  characteristics  o-f  the 
interleaved  FWT  (identical  components  at  successive  Cn]  ,  -for 
example).   Promising  matches  between  the  interleaved 
parameters  and  single  PRI  parameters  would  result  in  a 
decision  or  choice  o-f  which  single  PRI  trains  make  up  the 
interleaved  pulse  train. 

Scrutiny  o-f  the  PSD  plots  generates  similar  conclusions. 
Figure  4.16  shows  the  interleaved  PSD  plot.   Components  o-f 
relative  magnitude  o-f  9.3  or  greater  are  8,  13,  16,  19,  22, 
24,  25,  26,  and  32. 

Now  note  the  PSD's  o-f  TOA  TAG  5  and  8  (Figures  4.17  and 
D.29.    TOA  TAG  5  contributes  components  at  sequencies  o-f  13, 
22,  25,  and  26.  TOA  TAG  8  contributes  components  at  8,  16,  24, 
and  32. 

Be-fore  moving  on  to  another  interleaved  case,  a  discussion 
o-f  the  magnitudes  o-f  the  interleaved  components  (both  FWT  and 
PSD)  is  in  order.    A  de-finite  relationship  is  not  apparent 
between  components  o-f  individual  plots  and  interleaved  plots. 

The  author  can  say  that  one's  intuition  about  components 
o-f  interleaved  plots  generally  holds  true.   For  instance,  a 
large  positive  FWT  component  of  one  TAG  and  a  negative  FWT 
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Figure.  4.16.   PSD  :  TOA  TAG  5  AND  8  INTERLEAVED 
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Figure.    4,17.       PSD    :    TOA  TAG    5 
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component  o-f  another  TAG  at  the  same  Cn3  will  usually  generate 
a  smaller  component  in  a  FWT  component  of  the  interleaved 
TAG's.   The  result  is  not  a  simple  addition,  but  the  actual 
magnitude  o-f  the  interleaved  component  is  usually  smaller  than 
the  added  values  o-f  the  single  components.   No  -fast  and  easy 
rul e  appl ies. 

In  a  matching  algorithm,  the  magnitude  would  have  only  a 
comparative  bene-fit.   The  actual  promise  o-f  the  algorithm  lies 
in  matching  component  locations  (either  Cn3  or  sequency)  to  a 
data  bank  o-f  single  PRI  component  locations,  with  magnitude 
and  other  special  characteristics  as  a  decision  making  aid. 

One  more  interleaved  case  is  examined.   This  one  is 
different  because  there  are  no  coincident  pulses.   The  shift 
in  TOA  TAG  5  causes  it  to  interleave  perfectly  with  TOA  TAG 
19. 

TOA  TAG  5  SHF  1  (Figure  4.2  and  4.6)  DATA  line: 

DATA   8,1, 8, 9, 8, 9, 1,9, 9, 9, 0,1, 8, 8, 9, 9, 1,9, 9, 8, 9, 1,9, 9, 9, d, 1,8, 
8,8,8,1,8,8,8,8,1,8,8,8,8,1,8,8,8,8,1,8,8,8,8,1,8,8,9, 
9,1,9,9,9,9,1,9,9  (4.5) 

TOA  TAG  1?  (Figure  4.18  and  4.1?)  DATA  line: 

DATA   1,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,9,1,9,9,9,9,9,9,9,9, 
9,9,9,8,9,8,9,9,9,0,1,9,8,9,8,8,9,9,8,8,8,8,8,9,9,8,8, 
9,9,1,9,9,9,9,9,9  (4.6) 

And  the  interleaved  DATA  line: 

DATA   1,1,9,9,9,9,1,9,9,9,9,1,9,9,9,9,1,9,9,1,9,1,9,9,9,9,1,9, 
9,9,9,  1,9,9,9,9,  1,9,  1,9,9,  1,9,9,9,9,  1,9,9,9,9,  1,9,9,9, 
9,1,1,9,9,9,1,9,9  (4.5) 

Figure  4.29  and  4.21  show  the  FWT  and  PSD  of  the 

interleaved  TAG'S.   Interleaved  components  of  both  single  PRI 
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Figure.    4.18.       FWT    :    TOA  TAG    1? 
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Figure.    4.19.       PSD    :    TOA   TAG    1? 


00 

Si 

"  ™ 

J  00 

1  CM 

— i 

J  CD 

1  CM 

3^ 

1  CM 


4  CM 

1  CM 


*3 
CM 


J  00 

j  **  Hi 

^  a 

CM  CO 


Si 


i  oo 


CD 

CM 

SI 


SI 


159 


CO 


UJ 


< 
a: 

CO 

— 1 
< 

CD 

h- 

+ 

CO 

< 

IT) 

LL. 

in 


■q 


t 


■q 


si 


-L 


in 

1 


CD 

CD 

in 

CD 

SI 

c 

CM 
CO 

X 

LU 

a 

CM 

CO 

-J 
< 

CD 


CD 


SI 


€69992 '-I)    3VW  13d 


Figure.  4.29.   FWT  :  TOA  TAG  5  AND  1?  INTERLEAVED 


151 


-J 


CM 
00 

SI 
00 

GO 
CM 

CD 
CM 

CM 

CM 
CM 


I — 
CO 

2: 

LU 

a 

a 

LU 
> 

QC 

u 

21 
»— < 

0) 

+ 

a 

UO 

Q. 

si 

CM 


00 


>- 
u 
2: 

LU 

ID 

a 

LU 
4  CM    CO 


J  CO 

]      T-t 


S3 


GO 


CO 


—i  CM 


SI 


V*  O 


00 


co      in 


CO        CM 


Si 


(9S0Z0"=I)    3VW  138 


Figure.    4.21.       PSD    :    TOA   TAG    5   AND    19    INTERLEAVED 
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TAGS  are  present,  and  the  same  process  as  before  is  used  to 
identi-fy  them.   A  new  discovery  o-f  this  examination  reveals 
that  some  feature  PSD  components  o-f  TOA  TAG  5  SHF  1  are  passed 
directly  to  the  interleaved  PSD.   Note  the  magnitude  of 
sequencies   11,  16,  and  21  of  the  interleaved  PSD  is  the  same 
as  the  PSD  of  TOA  TAG  5  SHF  1,  and  that  the  magnitudes  of  TOA 
TAG  1?  PSD  at  these  sequencies  is  zero.   This  is  the  first 
evidence  of  two  single  train  representations  generating  an 
interleaved  components  magnitude  by  merely  adding  the  simple 
PRI  components.   A  more  closer  study  could  be  done  to  see  if 
it  is  merely  coincidence.   The  real  question  is,"  Did  the  fact 
that  there  are  no  coincident  pulses  in  the  interleaved 
representation  allow  me  to  add  certain  single  PRI  components 
to  form  an  interleaved  component?' 

A  partial  answer  to  the  question  probably  is  evident  since 
only  certain  components  of  the  two  single  PRI  PSD  plots  have 
this  feature  of  simple  addition. 

E.   CORRELATION  AND  THE  PSD  COEFFICIENTS 

Throughout  this  chapter  the  emphasis  has  been  placed  on 
the  location  and  relative  magnitudes  of  components  located  in 
the  Nal sh  index  domain  or  the  sequency  domain.   Most  of  the 
plots  show  the  magnitude  of  the  components  relative  to  the 
maximum  component  magnitude  of  that  particular  plot.   It  was 
stated  several  times  that  the  actual  magnitudes  of  the 
individual  PRI  (TOA  TAG)  components  didn't  seem  to  have  a 
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recognizable  relationship  with  the  magnitudes  of  the 
components  o-f  the  interleaved  plots. 

It  was  thought  that  the  magnitudes  might  be  related,  and  a 
correlation  o-f  the  interleaved  component  magnitudes  with  those 
o-f  the  individual  PRI  magnitudes  might  show  their  presence. 
The  PSD  coe-f-f icients  were  chosen  because  they  are  de-fined  -from 
the  FWT  coe-f-f  icients,  and  exhibit  the  same  characteristics  and 
information  about  the  trans-form  in  terms  o-f  symmetry  and 
relative  magnitude  components. 

Without  any  mathematical  -fan-fare,  the  correlation 

coefficient  P  is  presented. 

P  =    E    i    <X-X/>/tfx  •  <Y-Y/)/<*y  >       <4.6) 

where 

X,Y  are  random  variables 

X',  Y/    are  the  mean  of  the  variables 

<?x  y  are  *^e  variances  of  the  variables 

The  correlation  coefficient  expresses  the  degree  to  which 
two  random  variables  are  correlated  without  regard  to  the 
magnitude  of  either  one.   C38:  p. 83]    It  is  a  normalized 
quantity  and  will  be  a  value  between  -1  and  +1. 

In  this  application,  X  are  the  PSD  coefficients  of  the 
interleaved  situation  and  Y  are  the  PSD  coefficients  of  the 
individual  TAG  situation. 

A  correlation  BASIC  language  routine  was  written 
specifically  for  the  32  PSD  components  of  the  interleaved 
and  individual  TOA  TAGS.   This  program  was  adapted  from  a 
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correlation  routine  in  the  IMSL  library  contained  in  the  Naval 
Postgraduate  School  IBM  3833  computer.   It  was  not  adapted  -for 
general  application  but  speci-f ical  1  y  -for  this  situation.   For 
this  reason  it  is  not  provided  in  Appendix  A. 

The  correlation  routine  merely  reads  in  the  component 
values  o-f  the  individual  PRI  TAG  representations  and  the 
interleaved  representation  generated  by  the  FWT/PSD/Group 
Spectrum  program  provided  in  Appendix  A.   It  correlated  the 
two,  based  upon  Equation  (4.6) ,  and  outputed  a  single  number, 
the  correlation  coefficient. 

The  PSD  coefficients  of  the  interleaved  situation  of  TOA 
TAG  5  and  TOA  TAG  1?  were  used  (Figure  E.5  in  Appendix  E) . 
This  interleaved  situation  is  different  from  that  of  page  148 
because  TOA  TAG  5  wasn't  shifted  before  interleaving,  causing 
the  first  pulse  of  TOA  TAG  5  and  19  to  be  coincident.  With 
these  coefficients,  the  values  of  the  individual  TOA  TAG 
components  from  3  to  64  were  then  correlated  with  the 
interleaved  values  using  the  correlation  routine. 

Table  4.2  gives  the  TOA  TAG  number,  and  the  correlation 
coefficient  generated. 

Examination  of  the  correlation  coefficients  shows  strong 
correlation  with  TOA  TAG  5.   Great!   It  is  desired  to  have  a 
high  correlation  value  with  TOA  TAG's  5  and  19,  since  the 
interleaved  representation  was  composed  of  these  two  TAG 
representations.   The  TOA  TAG  5  correlation  coefficients  fits 
the  need  nicely. 
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Table  4.2.   Correlation  o-f  Interleaved  PSD  coe-f-f  icien  ts  with 

individual  TOA  TAGS 


TOA  TAG  number  Correlation  Coe-f -f icien t 


3  8.665 

4  8 . 234 

5  8.963 

6  8.582 

7  8.385 

8  8.838 

9  8.166 

18  8.789 

11  8.299 

12  8.183 

13  8.274 

14  8.191 

15  8.365 

16  -8.865 

19  8.148 
22  8.338 
25  -8.828 
32  -8.855 
48  8.186 
51  8.859 
64  -8.722 


On  the  down  side,  the  interleaved  PSD  coe-f  -fie  ients  do  not 
show  a  high  value  o-f  correlation  with  TOA  TAG  19.  The  initial 
success  o-f  TOA  TAG  5  certainly  indicates  that  this  is  an  area 
■for    -further    investigation,    however. 

The    author    suggests    that    the   correlation    procedure   could 
be   used    in    an    e-f-fective    deinterl  eaving    algorithm    i-f    it   was 
applied    to    the    sequential    method   o-f    extracting   PRI's.    Once    the 
initial    correlation    routine    is   carried   out    and   a   high   value   o-f 
correlation    coe-f-f icient    is    identi-fied    and   associated  with    an 
individual    TOA  TAG   string,     then    the    algorithm   could   return    to 
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the  interleaved  DATA  stream  representation  and  subtract  (set 
equal  to  zero)  those  TOA  TAG'S  associated  with  the  determined 
PRI  . 

Now  repeat  the  entire  process.   Take  the  FWT  and  PSD  o-f 
the  remaining  TOA  TAG  representation,  correlate  the  PSD 
coeff icients  with  the  data  bank  o-f  individual  TOA  TAG  PSD 
coeff icients,  and  see  i-f  a  new  TOA  TAG  has  been  identi-fied 
with  another  high  value  o-f  correlation  coe-f  f  i  cient .   I-f  it 
has,  then  strip  these  TOA  TAG'S  -from  the  original  string. 

The  process  could  be  repeated  until  calculation  o-f  the  PSD 
yields  no  -further  use-ful  coe-f-f  icients,  or  until  no  correlation 
IS  -found.   The  stripped  TAG's  represent  the  individual  pulse 
trains  that  make  up  the  interleaved  pulse  string. 

This  suggested  process  is  similar  to  the  way  a 
preprocessor  sorts  the  incoming  data  stream.   A-fter 
identifying  particular  pulses  in  the  stream,  the  preprocessor 
strips  these  pulses  -from  the  stream,  simplifying  the  -follow-on 
process,  which  continues  to  identify  <dei n ter 1 eave)  individual 
pulse  trains  of  the  remaining  stream.   Each  identified  train 
is  then  removed  from  the  incoming  data  from  the  receiver. 

A  problem  that  plagues  current  deinter 1 eaving  techniques, 
that  of  coincident  or  missing  pulses,  would  probably  plague 
this  correlation  deinter 1 eaving  procedure. 

The  TOA  TAG  representation  does  not  indicate  coincident 
pulses  separately.   Coincident  pulses  are  represented  by  a 
single  TOA  TAG.   Only  the  receiver  knows  which  one  it  detected 
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and  represented  as  a  received  pulse.   There-fore,  the  stripping 
of  this  TAG  generated  when  two  pulses  had  arrived  at  the 
receiver  at  the  same  time  inadvertently  strips  a  TOA  TAG  -from 
the  remaining  pulse  representation. 

I-f  the  PSD  coeff icients  of  the  remaining  string  with  the 
inadvertently  stripped  TAG  were  now  calculated,  one  would  find 
they  are  different  than  those  calculated  from  the  remaining 
string  that  included  the  inadvertently  stripped  TAG. 

These  PSD  coefficients,  those  calculated  from  the  missing 
TAG  string,  cannot  correlate  higher  with  the  individual  TOA 
TAG  strings  (used  by  the  correlation  routine)  than  the 
coe-ff  icients  o-f  the  complete  remaining  string. 

The  individual  TAG  PSD  coefficients  were  calculated  from 
TOA  TAG  strings  with  no  missing  pulses.   So,  coefficients  from 
strings  with  missing  pulses  will  not  correlate  as  well  as 
coe-ff  icients  calculated  from  complete  strings. 

A  possible  solution  could  be  the  use  of  a  large  number  of 
TOA  TAG's  in  the  interleaved  representation,  say  1824  or  2848. 
In  other  words,  expand  the  "window"  or  aperture  on  the  data 
stream.  This  would  allow  the  interleaved  representation  to  be 
more  realistic  and  comparable  to  the  receiver's  own  TOA 
resolution.   More  coefficients  might  help  the  correlation 
procedure  and  lessen  the  effect  of  missing  and  coincident 
pu 1 ses . 
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g*   CONCLUSIONS  AND  RECOMMENDATIONS 

Be-fore  discussing  the  conclusions  of  the  research  e-f-fort 
and  making  recommendations  -for  -follow-on  work,  a  review  o-f  the 
investigation  steps  that  were  taken  by  the  author  and 
translated  into  the  substance  o-f  this  thesis  will  be 
beneficial  . 

For  background,  the  importance  o-f  the  ESM  e-f-fort  and  its 
place  in  Electronic  War-fare  was  studied.   A  typical  ESM 
system  was  examined,  and  the  problem  of  deinterl eaving  the 
pulse  stream  output  from  an  ESM  receiver  was  considered,  along 
with  current  techniques  for  deinter 1 eaving  already  in  use. 

Having  no  previous  idea  of  even  the  existence  of  the  Walsh 
functions  set  the  stage  for  the  study  of  these  unique  and 
interesting  functions.   As  the  study  proceeded  through  the 
definition  of  the  Walsh  Transform,  the  constant  comparisons  of 
the  Walsh  Transform  properties  with  the  familiar  properties  of 
the  Fourier  Transform  brought  review  of  this  area. 

Reference  23  provided  the  FWT  program  that  became  the  tool 
of  the  investigation.   Adaptation  of  the  program  to  the  BASIC 
language  and  a  gathering  of  evidence  that  supported  belief  in 
its  output  took  a  considerable  time.   Searches  of  the 
literature  for  examples  of  FWT's  of  familiar  waveshapes 
finally  provided  the  proof  of  the  program. 

With  the  guidance  of  the  thesis  advisor,  the 
representation  of  the  receiver  data  stream  with  time  of 
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arrival  tags  was  accepted,  and  the  FWT  program  was  applied  to 
many  individual  tag  strings.   The  purpose  was  to  determine  if 
a  TOA  TAG  string  representing  a  particular  PRI  pulse  train 
could  be  recognized  by  certain  features  of  the  Nal sh  Transform 
or  the  Power  Spectral  Density.   The  FWT,  PSD,  and  Group 
Spectrum  of  interleaved  TOA  TAGS  were  then  closely  examined 
for  PRI  recognition  features.   Finally,  a  correlation 
algorithm  was  used  to  determine  the  degree  of  correlation 
between  interleaved  PSD  coeffcients  and  single  pulse  train  PSD 
coefficients.   An  examination  of  the  correlation  coefficients 
was  done  to  see  if  a  PRI  could  be  determined  by  noting  the 
magnitude  of  the  coefficient. 

A.   CONCLUSIONS 

One  result  of  the  effort  was  the  demonstration  of  the 
ease  of  computation  of  the  Walsh  Transform  and  Power  Spectral 
Density  of  a  number  series.   The  algorithm  is  fast  and 
efficient,  and  easily  adapted  to  a  BASIC  language 
implementation.   This  allows  the  computation  of  the  transform 
and  calculation  of  the  PSD  coeffcients  on  a  microcomputer, 
easily  used  and  readily  available  in  the  work  space. 

When  applied  to  a  TOA  TAG  representation  of  a  single  PRI 

pulse  train,  the  transform  shows  features  that  can  be  used  to 

recognize  this  pulse  train  when  compared  to  the  transforms  of 

other  PRI  pulse  trains.   Symmetry  of  the  coefficient 

distribution,  and  the  relative  magnitude  of  certain  Walsh 

index  or  sequency  components  are  the  features  that  allow  this 

identification . 
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It   was   shown    in    Chapter    4    that    each    single    PRI    pulse 
train's      Nalsh    trans-form   has    unique    and    distinctive    components 
that    distinguish    it    -from   other    transforms   o-f    single   PRI 
representations.      No    single    feature,    however,    such    as    a 
relative  maximum   sequency    component,    can    be   connected   with    the 
PRI    of    the    pulse    train.       Each    particular    PRI    transform    simply 
has    individual    features    that    identify    it,    and   only    certain 
characteristics,    such    as   symmetry,    are    common    between    groups 
of    individual    pulse    train    transforms. 

The  Walsh    transform    and   PSD's    of    interleaved  TOA   TAG 
representations   have   features    that    can    be    identified    as 
belonging    to    the    transforms    or    PSD's   of    component    individual 
pulse    trains.       However,    these    features    are   not    readily 
apparent    and   one   could    not    identify   or    associate    them  with    a 
particular    single   PRI    transform  without    prior    knowledge   of 
single   PRI    transform    features.       Like    the    transforms   of    single 
PRI    pulse    trains,     the    interleaved    transform   or    PSD    does   not 
have    any   features    that    can    identify    component    PRI's    from    the 
transform   or    PSD    alone. 

This    lack    of    component    PRI    features    in    the    transform    or 
PSD    effectively   prevents    its   use    in    a   practical    deinter 1 eaving 
algorithm.      At    the   outset    of    the   research    effort,     it   was    hoped 
that    the    transform    of    a   representation    of    the   pulse   stream 
from   a   ESM   receiver   would    identify    the    individual    PRI's    of    the 
component    pulse    trains.      Neither    the    transform   or    the   PSD   does 
this,    although    there    are   common    features    between    the 
interleaved    transform/PSD    and   component    transform/PSD. 
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The  application  o-f  a  correlation  routine  between  the 
the  interleaved  PSD  coeff icien ts  and  the  individual  TOA 
TAG  PSD  coefficients  was  seen  to  be  effective  in  identifying 
one  o-f  the  component  PRI's.   It  did  not  identify  all 
component  PRI's,  and  any  algorithm  that  uses  this 
approach  will  have  to  go  through  additional  steps  to  do  so. 
These  required  additional  steps  quickly  negate  the  speed  o-f 
computation  advantage  that  the  Walsh  trans-form  exhibits,  and 
the  approach  becomes  quite  similar  to  current  sequential  PRI 
extraction  techniques. 

In  summary,  the  transform  o-f  the  interleaved  TOA  TAG 
representation  -failed  to  produce  a  simple  PRI  recognition 
•feature.   Recognition  -features  are  present,  howler ,  in  the 
interleaved  trans-form  and  PSD  that  would  allow  the 
identification  of  component  pulse  trains  in  an  interleaved 
representation  if  a  data  banK  of  individual  PRI  transform/PSD 
features  was  available. 

B.   RECOMMENDATIONS 

It  is  recommended  that  an  effort  at  writing  a  "matching 
feature"  deinter 1 eaving  program  be  made.   The  program  would 
use  a  data  bank  of  parameters  that  describe  the  interesting 
properties  and  features  of  the  Walsh  transform  and/or  PSD  of  a 
TOA  TAG  representation  of  a  single  PRI  pulse  train. 

The  data  bank  would  be  the  measure  of  attaining  a 
successful  deinter 1 eaving  program.   Properties  such  as  even 
symmetry,  double  components  of  power  of  two  interval 
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representations,  and  distinctive  index  or  sequency  components 
of  the  odd  numbered  TOA  TAG  representations  would  be 
included. 

The  algorithm  would  scan  the  FWT/PSD  of  the  interleaved 
TOA  TAG  representation  for  recognizable  features,  then  match 
these  features  to  the  data  bank.   Appropriate  tolerances  could 
be  included  to  be  used  in  decisions  about  whether  a  single  PRI 
is  a  component  of  the  interleaved  train. 

Additional  study  could  possibly  yield  the  optimal  features 
to  include  in  the  bank.   For  certain,  a  closer  look  at  single 
PRI  representations  must  be  done  to  determine  which  components 
are  shared  among  them.   The  author  firmly  believes  that  there 
will  always  be  a  FWT/PSD  feature  of  a  single  PRI 
representation,  whether  a  particular  index  or  sequency 
component,  or  a  symmetrical  characteristic,  that  can  be  used 
to  distinguish  it  from  other  single  PRI  representations. 

Closer  study  of  whether  there  is  a  relationship  between 
interleaved  and  individual  index  or  sequency  magnitudes  should 
be  completed.   It  is  a  feature  that  could  possibly  be  used  in 
a  matching  algorithm,  and  it  certainly  can  be  used  in  direct 
compuation  methods,  such  as  a  least  mean  squares  or 
correlation  routine  that  might  generate  an  individual  PRI 
indication . 

The  quite  tentative  but  interesting  result  found  in  the 
last  section  of  Chapter  4  should  be  explored  completely.   All 
coefficients  should  be  involved,  including  the  Group  Spectrum 
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coe-f-f icients.   In  this  thesis,  the  Group  Spectrum  coefficients 
have  mainly  been  used  to  show  the  time  invariance  nature  of 
their  character,  but  they  might  offer  some  interesting  direct 
computation  features. 
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APPENDIX  A 
FAST  WALSH  TRANSFORM  PROGRAMS 

The  following  pages  contains  listings  o-f  the 
Walsh  Trans-form  programs  used  in  this  thesis.   The  -first 
program  is  written  -for  an  IBM  Personal  Computer,  and  the 
second  was  written  -for  the  HP-85 . 

The  FWT  subroutine  was  adapted  from  a  FORTRAN  program  in 
■Walsh  Functions  and  their  Applications,"  by  K.G.  Beauchamp. 

The  PSD  and  Group  Spectrum  subroutines  were  written  by  the 
author  . 

Both  programs  were  run  on  the  respective  computer  BASIC 
interpreter  < uncompi 1 ed) .   A  compiled  version  would  probably 
run  faster,  although  a  128  number  FWT  ran  on  the  IBM  PC  in 
about  4  seconds  using  only  the  interpreter. 
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18  '   This  program  computes  the  Fast  Walsh  Tran-form  ,  the 

Power  Spectral  Density,  and  the  Group  Spectrum 
20  '   Coe-f-f icients  o-f  a  number  series. 

38  '   The  FWT  subroutine  is  a  BASIC  adaptation  o-f  a  FORTRAN 
program  listed  48  '  in  "  Walsh  Functions  and  Their 
Applications,"  by  K.G.  Beauchamp. 
'   The  PSD  and  Group  Spectrum  subroutines  were  written  by 

the  author . 
'   This  program  runs  on  a  IBM  Personal  Computer.   The 

printer  used  was  a  NEC  8023. 
DIM  X< 100)  ,Y<64)  ,P<64)  ,G<  11)  ,F(64)     'dimension  arrays 
LPRINT  CHR*C27)  ;'LU  ;"008"      'set  le-ft  margin 
N  =  64  '  number  o-f  input  numbers 

FOR  I  =  1  TO  N       '  read  in  data  

'READ  X<I)         '-from  the  DATA  line,  or 

'calculate  your  own 


45 

50 

60 
65 
70 
80 
90 

100 

110 
120 
130 
140 
150 
160 
170 
188 
190 

200 

230 
235 
240 
250 
260 

270 
280 
290 

300 

310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 


o-f  input  series 

coe-f-fs  by  N  and  downsh-ft 

contains  the  FWT  coe-f-fs 


XU)  =  SIN( 1XI/N*2£3. 1416) 
NEXT  'I 

GOSUB  310  'calculate  FWT 
FOR  B  -  1  TO  N  'divide  all 
F(B-l)  =  X<B)/N  '  this  array 
NEXT  'B 

'calculate  power  spectrum  coe-f -f  ici  ents 
'  calculate  group  spectrum 
'output  the  FWT  coe-f-fs 
'output  the  PSD  coe-f-fs 
'output  the  group  coe-f-fs 


GOSUB 
GOSUB 
GOSUB 
GOSUB 
GOSUB 
END 


600 

720 
840 
950 
1047 


SUBROUTINE  FHT<N,X,Y) 


in  array 


This  routine  per-forms  a  FWT  of  an  input  series 
X.   The  array  Y  is  used  -for  working  space. 
The  dimensions  o-f  X  and  Y  must  be  a  power  o-f  2. 
The  results  o-f  the  FWT  are  in  sequency  order,  positive 
phasing,  and  in  array  X.  The  subroutine  uses  a  Hadamard 
trans-form. 


N2  =  N/2 
M  -  6 
FOR  L  =  1  TO  M 
NY  =  0 
NZ  -  2-ML-l) 
NZI  -  2  *    NZ 
NZN  -  N/NZI 
FOR  I  *  1  TO  NZN 
NX  =  NY  +  1 
NY  =  NY  +  NZ 
JS  =  U-l)  X  NZI 
JD  ■  JS  +  NZI  +  1 
FOR  J  =  NX  TO  NY 
JS  =  JS  +  1 


M  equals  the  base  2  logarithm  o-f  N 
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458        J2  =  J  +  N2 

468        Y<JS)  -  X<J)  +  X<J2> 

478        JD  =  JD  -  1 

488        Y<JD>  ■  X(J>-  X<J2> 

498      NEXT  'J 

588    NEXT  'I 

518   FOR  B  -  1  TO  N 

528   X<B)  =  Y<B> 

538   NEXT  'B 

548  NEXT  'L 

558  RETURN 

555  '  SUBROUTINE  PSD<N) 

568  ' 

578  ' 

588  '   This  routine  calculates  the  Walsh  Power  Spectral 

Density  Coef f icien ts. 
598  ' 

688  P<8)  =  <F<8))A2  '-first  PSD  coe-f f ici en t 

618  K=l 

628  FOR  I  -  1  TO  <N/2)-l     'coefficients  1  to  N/2-1 
638      P<I)=  (F<K>>*2  +  <F<K+1>)~2 
648      K  -  K  ♦  2 
658  NEXT  'I 

668  P<N/2>=<F<N-1>>  A2      'last  PSD  coefficient 
678  RETURN 
688  ' 

698  '  SUBROUTINE  GROUP  SPECTRUM(M) 

788  ' 
718  '  This  subroutine  calculates  the  group  spectrum 

coefficients. 
728  D=l 

738  G<8>  -  P<8)A2 
748  FOR  B  =  8  TO  M-l 

758  FOR  C  =<2^B)  TO  N/2-1   STEP  <  2*2'KB> 
768  G<D)  =  G<D)  +  P<C)A2 
778  NEXT  'C 
788  D  ■  D+l 
798  NEXT  'B 
888  G(M>  -  P<N/2>*2 
818  RETURN 
828  ' 
838  '   This  subroutine  outputs  the  Walsh  Transform 

Coefficients. 
848  ' 

841  FOR  C  =  1  TO  5:  LPRINT:  NEXT       '  1  inch  top  margin 

842  LPRINT"  FWT  COEFFICIENTS  FOR  1  CYCLE  SINE 
FUNCTION" 

843  LPRINT:   LPRINT    'skip  two  lines 

858  LPRINT  "Walsh  Function"  TABC28)  "Coeff."  TAB<35) "Wal sh 
Function"  TAB<55)  "Coeff." 
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868  LPRINT  STRING*< 15, 45)  TAB<28)  STRING*< 6 ,45)  TAB<35) 

STRING*< 15,45)  TAB<55)  STRING*< 6,45) 
878  FOR  B  ■  1  TO  32      'output  the  coe-f-fs 


B-l;",t)"  TAB< 16) :   LPRINT  USING 
F(B-l);:  LPRINT"  ■     TAB<35) 

B+31;",t)"  TAB<52)  :   LPRINT  USING 
F<B+31) 


888  LPRINT  "WAL< " 

898  LPRINT  "WAL< " 

988  NEXT  'B 

918  FOR  C  =  1  TO  16  :  LPRINT:  NEXT     'skip  16  lines 

911  LPRINT0  xx"    'print  pg.  no 

912  LPRINT  CHR*< 12)   '-form  -Feed 
928  RETURN 

938  ' 

948  '        Output  the  PSD  coe-f  -f  i  ci  en  ts . 

958  ' 

951  FOR  C  =  1  TO  5:  LPRINT:  NEXT 

978  LPRINT  "  PSD  COEFFICIENTS  FOR  A  1  CYCLE  SINE 

FUNCTION" 
988  LPRINT  :  LPRINT 
998  FOR  B  =  8  TO  N/2 

1888  LPRINT  "P<";  B  ;")  =  "  ;:  LPRINT  USING  " tt* .#*#**##"; P< B) 
1828  NEXT  'B 
1825  LPRINT:  LPRINT 
1848  RETURN 

1844  ' 

1845  '   Output  the  group  coe-f-f icients 

1846  ' 

1847  LPRINT"     GROUP  SPECTRUM  COEFFICIENTS  OF  A  1  CYCLE  SINE 

FUNCTION" 

1848  LPRINT:  LPRINT 
1858  FOR  B  =  8  TO  M 

1868  LPRINT  "G(";B;")=  "; -.LPRINT  USING  "#*.  *#*##*#"  ;  G<  B) 
1878  NEXT  'B 

1871  FOR  C  =  1  TO  7:  LPRINT:  NEXT 

1872  LPRINT  "  xx" 

1873  LPRINT  CHR*<  12)     '-form  feed 
1888  RETURN 
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18  !   This  program  computes  the  Fast  Walsh  Tranform  ,  the 

Power  Spectral  Density,  and  the  Group  Spectrum 
28  !   Coefficients  of  a  number  series. 
38  !   The  FWT  subroutine  is  a  BASIC  adaptation  of  a  FORTRAN 

program  listed  48  !  in  "  Walsh  Functions  and  Their 

Applications,"  by  K.6.  Beauchamp. 
45  !   The  PSD  and  Group  Spectrum  subroutines  were  written  by 

the  author. 
58  !   This  program  runs  on  HP  85  Computer,  and  writes  the 

coefficients  to  data  files  named  "FWT" , "PSD" , and  "GRP" 
48  DIM  X< 188) ,Y<64) ,P<64) ,G< 11) ,F<64)     Idimension  arrays 
65  PRINTER  IS  2     !  select  printer  for  output  of  print 
70  N  =  64  !  number  of  input  numbers 

88  FOR  I  ■  1  TO  N       !  read  in  data  

98  '.READ  X<I)         !from  the  DATA  line,  or 

188  X(I)  =  SIN< 1XI/N*2*3.  1416)    Icalculate  your  own 

118  NEXT  I 

128  GOSUB  318       Icalculate  FWT  of  input  series 

138  FOR  B  =  1  TO  N     ! divide  all  coeffs  by  N 

148  X<B)  -  X<B)/N   !  this  array  contains  the  FWT  coeffs 

158  NEXT  B 

155  BEEP  Q    CLEAR 

156  DISP  "FWT  CAL  OVER" 

168  GOSUB  688     '.calculate  power  spectrum  coefficients 

165  DISP  "PSD  CAL  OVER" 

178  GOSUB  728     !  calculate  group  spectrum 

175  DISP  "GRP  COEFF  CAL  OVER" 

188  GOSUB  348   ! output  the  FWT  coeffs 

198  GOSUB  958   ! output  the  PSD  coeffs 

288  GOSUB  1847   ! output  the  group  coeffs 

285  CREATE  "FWT", 38         !  creates  file 

286  ASSIGN  *1  TO  "FWT'    !  opens  file 

287  FOR  B  =  1  TO  64 

288  PRINT*  1;  X<B) 

289  NEXT  B 

218  ASSIGN*  1  TO  X       ! c 1 oses  file 

211  CREATE  "PSD" ,  28 

212  ASSIGN  *1  TO  "PSD' 

213  FOR  B  =  8  TO  N/2 

214  PRINT*  1;  P<B) 

215  NEXT  B 

216  ASSIGN*  1  TO  * 

217  CREATE  "GRP* ,  28 

218  ASSIGN  *1  TO  "GRP' 

219  FOR  B  -  0  TO  M 
228  PRINT*  1;  G(B) 

221  NEXT  B 

222  ASSIGN*  1  TO  * 

225  BEEP  @  DISP  "PROGRAM  OVER" 
238  END 
231  ! 
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235 
248 
258 
268 

278 
288 
298 

388 
318 
328 
338 
348 
358 
368 
378 
388 
398 
488 
418 
428 
438 
448 
458 
468 
478 
438 
498 
588 
518 
528 
538 
548 
558 
555 
568 
578 
588 

598 
688 
618 
628 
638 
648 
658 
668 
678 
671 
672 


SUBROUTINE  FHT<N,X,Y) 

This  routine  performs  a  FWT  of  an  input  series  in  array 
X.   The  array  Y  is  used  -for  working  space. 
The  dimensions  of  X  and  Y  must  be  a  power  of  2. 
The  results  of  the  FWT  are  in  sequency  order,  positive 
phasing,  and  in  array  X.  The  subroutine  uses  a  Hadamard 
transform. 


N2  =  N/2 
M  =  6 
FOR  L  =  1  TO  M 
Nl  =  8 
N3  =  2"<L-D 
N4  =  2  *  N3 
N5  =  N/N4 
FOR  I  -  1  TO  N5 
N6  =  Nl  +  1 
Nl  =  Nl  +  N3 
Jl  ■  <I-1)  *  N4 
J3  -  Jl  +  N4  +  1 
FOR  J  =  N6  TO  Nl 
Jl  -  Jl  +  1 
J2  =  J  +  N2 
Y<J1)  ■  X(J)  +  X<J2) 
J3  =  J3  -  1 
Y<J3)  =  X<J>-  X<J2) 
NEXT  J 
NEXT  I 
FOR  B  =  1  TO  N 
X<B>  =  Y<B> 
NEXT  B 
NEXT  L 
RETURN 


M  equals  the  base  2  logarithm  of  N 


SUBROUTINE  PSD<N) 


This  routine  calculates  the  Walsh  Power  Spectral 

Density  Coefficients, 
i 

P<8)  -  <F<8>)*2  Ifirst  PSD  coefficient 

K=l 

FOR  I  =  1  TO  <N/2)-l     !  coefficients  1  to  N/2-1 

P<I)  =  <F(K)>A2  +  <F<K+1)>*2 

K  -  K  +  2 
NEXT  ! I 

P<N/2)=<F<N-1>)  A2      ".last  PSD  coefficient 
RETURN 
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688 
690 
788 
718 


SUBROUTINE  GROUP  SPECTRUM<M) 


This  subroutine  calculates  the  group  spectrum 

coefficients. 
728  D=l 

738  G<8)  ■  P<8)A2 
748  FOR  B  =  8  TO  M-l 

758  FOR  C  =<2~B)  TO  N/2-1   STEP  <2*2AB) 
768  6(D)  =  G<D)  +  P<C)*2 
778  NEXT  !C 
788  D  -  D+l 
798  NEXT  B 
888  G<M)  -  P<N/2)A2 
818  RETURN 
828  ! 
838  !   This  subroutine  outputs  the  Walsh  Trans-form 

Coefficients. 
848  ! 

858  PRINT  "Walsh  Function" ;TAB< 28) ; "Coeffi ci en t"  @  PRINT 
878  FOR  B  =  1  TO  64      ! output  the  coeffs 
888  PRINT  "WALC  ;  B- 1 ;  "  ,  t)  "  ;TAB<28)  ;  F(B) 
988  NEXT  B 
928  RETURN 
938  ! 

948  !        Output  the  PSD  coefficients. 
958  ! 
968  PRINT 

978  PRINT  "POWER  SPECTRUM  COEFFS"  6  PRINT 
998  FOR  B  =  8  TO  N/2 
1888  PRINT  "P<";  B  ;")="  ;  P<B) 
1828  NEXT  B 
1848  RETURN 

1844  ! 

1845  !   Output  the  group  coefficients 

1846  ! 

1848  PRINTS  PRINT 

1849  PRINT  "GROUP  SPECTRUM  COEFFICIENTS"  S  PRINT 
1858  FOR  B  =  8  TO  M 

1868  PRINT  "G<";B;  ")=  ";@PRINT  USING  ■»#.«**##**"; G<B) 

1878  NEXT  B 

1888  RETURN 

1885  !   This  is  a  typical  DATA  line  of  numbers  representing 

TOA  TAGS 
1898  !DATA  1,8,8,8,8,8,8,8,1,8,8,8,8,8,8,8,1,8,8,8,8,8,8,8,1, 
8,8,8,8,8,8,8,1    'etc.  for  N  =  64 
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APPENDIX  B 
ORTHOGONALITY 

The  concept  of  orthogonality  in  a  set  o-f  functions  is 
important  because  only  orthogonal  sets  o-f  functions  can  be 
made  to  represent  another  function  with  a  required  degree  of 
accuracy.   The  term  itself  often  brings  into  mind  the  word 
"perpendicular',  and  this  thought  can  be  a  visualization  of 
the  structure  of  the  members  of  an  orthogonal  set. 

Consider  that  we  have  a  function  set,  Sn(t) ,  where  n  = 

0,1,2,3,...    The  set  is  said  to  be  orthogonal  with  weight  K 

over  the  interval  8>t>T  if 

T  K   i  f  n=m 

K-Sn(  t)Sm<t)dt  =  <B.  1) 

8  8      n*«m 

with  n  and  m  being  integer  values.   If  the  constant  K  is  equal 

to  1  then  the  set  is  normalized  and  the  set  is  referred    to  as 

an  orthonormal  set. 

With  this  orthogonal  set,  we  may  now  represent  another 

function,  f<t),  defined  over  the  interval  (8,T),  as 

T 

8 


f(t)  =  ZCnSn<t)  <B.2> 


and  Cn  is  a  number  that  indicates  the  value  of  the  function 
Sn<t) .   Cn  can  be  chosen  so  as  to  mimimize  the  mean-square 
error  in  representing  f<t). 
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Cn  -  1/T  •   -f<t)Sn<t)dt  <B.3> 


T 
■f 
9 

This  orthogonal  -function  series  representation  reduces  the 

number  of  coefficients  needed  to  completely  represent  the 

signal.   C23:  pp. 1-3] 
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APPENDIX  C 
M0DUL0-2  ADDITION  AND  THE  GRAY  CODE 

Modulo-2  addition  is  an  important  mathematical  operation 
in  Walsh  Theory.  Its  operation  is  used  in  the  de-finition  o-f 
the  Paley  ordered  Walsh  -functions  with  Rademacher  -functions, 
and    also    in    the   product    o-f    two  Walsh    -functions,    namely, 

WAL<n,t)WAL<m,t>    =WAL<n«m,t>  <C.1> 

where   n$m    indicates  modulo-2    addition. 

Modulo-2    addition    is   really   binary    sums   without    the    ca^rry, 
and    obey    the   rules 

8+8=8,  8+1=1,  1+8=1,    and       1+1=8 

The  Gray  Code  is  a  binary  code  that  is  o-ften  used  in 
communications  because  the  codes  -for  successive  decimal  digits 
di-f-fers  by  only  1  bit.   It  is  not  a  weighted  code,  meaning 
that  the  decimal  value  o-f  a  coded  digit  cannot  be  computed  by 
a  simple  -formula. 

Table  C.l.   Gray  Code  -for  16  Digits. 
Decimal         Code  Decimal         Code 


8 

8888 

1 

888  1 

2 

88  11 

3 

8818 

4 

8118 

5 

8111 

6 

8181 

7 

8  188 

8 

1188 

9 

118  1 

18 

1111 

11 

1118 

12 

18  18 

13 

1811 

14 

1881 

15 

1888 
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APPENDIX  D 
FAST  WALSH  TRANSFORM  AND  POWER  SPECTRAL  DENSITY  PLOTS 

The  -following  pages  are  the  plots  o-f  the  FWT  and 
PSD  coe-f-f icient  o-f  the  simulated  TOA  TAG  strings  used  in  the 
thesis. 

The  plots  are  normalized  to  the  maximum  component  and  thus 
show  the  relative  value  o-f  each  component  to  this  maximum. 

The  value  o-f  the  maximum  component  is  indicated  on  the 
vertical  axis. 
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Figure   D.i.       FWT    :    TOA  TAG   2 
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Figure    D.2.       FWT    :    TOA   TAG    3. 
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Figure  D.3.   FWT  :  TOA  TAG  4 
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Figure    0.4.       FWT    :    TOA   TAG    5 
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Figure    D.5.       FWT    :    TOA   TAG    6 
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Figure    D.6.       FWT    :    TOA   TAG    7. 
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Figure    D.7.       FWT    :    TOA   TAG    8 
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Figure    D.8.       FWT    :    TOA   TAG    ? 
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Figure    D.9.       FWT    :    TOA   TAG    19 
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Figure    D.18.       FWT    :    TOA   TAG    11 
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Figure    D.U.       FWT    :    TOA   TAG    12 
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Figure    D.12.       FWT    :    TOA   TAG    13 
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Figure    D. 13.       FWT    :    TOA   TAG    14 


CO 

CD 

00 

la 

CM 

en 

X 

UJ 

cm 

X 

en 
< 

CD 


CD 


3 


138 


Li- 
en 

< 


CO 


in 


CD 

< 

< 

a 


CO 

< 


IT) 


S3 


-L 


in 

! 


ci8Z.0"=i)  dvw  lay 


Figure    D.14.       FWT    :    TOA   TAG    15 
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Figure    D.15.       FWT    :    TOA   TAG    16 
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Figure    D.16.       FWT    :    TOA   TAG    19. 
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Figure   D.17.       FWT    :    TOA  TAG   22. 


192 


LL- 
CO 

-z. 

<  in 

££   CM 

i— 

X   < 

CO    h- 

-J 

<  < 

2t  a 


CO 

< 

Li- 


~ I 


■r 


CD 

CD 

in 

CD 

V 

C\i 
00 

X 

Ui 

a 
*—* 

CO 

CD 


ao 


s 


m 


cs 


in 


(69*0 '-!)    DVW  13d 


Figure    D.13.       FWT    :    TOA   TAG    25 
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Figure    D.19.       FWT    :    TOA   TAG    32. 
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Figure  D.28.   FWT  :  TOA  TAG  49 
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Figure    D.21.       FWT    :    TOA   TAG    51. 
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Figure    D.22.       FWT    :    TOA   TAG    64. 
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Figure   D.23.       PSD    :    TOA  TAG   2 
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Figure    D.24.       PSD    :    TOA   TAG    3 
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Figure    D.25.       PSD    :    TOA   TAG    4. 
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Figure    D.26.       PSD    :    TOA   TAG    5. 
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Figure    D.27.       PSD    :    TOA   TAG    6 
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Figure    D.28.       PSD    :    TOA   TAG    7 
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Figure    D.29.       PSD    :    TOA   TAG    8. 
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Figure    D.38.       PSD    :    T0A   TAG    9 
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Figure    D.31.       PSD    :    TOA   TAG    19 
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Figure    D.32.       PSD    :    TOA  TAG    11. 
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Figure    D.33.       PSD    :    T0A   TAG    12. 
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Figure    D.34.       PSD    :    TOA   TAG    13 
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Figure    D.35.       PSD    :    TOA   TAG    14 
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Figure    D.36.       PSD    :    TOA   TAG    15 
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Figure  D.37.   PSD  :  TOA  TAG  16 
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Figure    D.38.       PSD    :    TOA   TAG    1?. 
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Figure    D.39.       PSD    :    TOA   TAG    22 


H  OJ 


SI 


*-•        O) 


a> 


CO 


in 


CO 


oj      »-* 


is 


214 


CM 
CO 


CO 

zz.  in 

Ui    CM 

a 

CD 

cj  < 

£L 
CO 


UJ 

a 


<: 
a 


SI 

00 

CM 

CO 
CM 

CM 

CM 
CM 

Si 
CM 

CD 
»-« 

CO 


i  SI 

i 

-4 

-I  CO 


C_3 


UJ 

cm  en 


4  CM 


SI 


•-*        O)        CD 


co      in 


CO       cm       ^ 


(W!0-=l)    3VW  "Bed 


Figure  D.48.   PSD  :  TOA  TAG  25 
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Figure    D.41.       PSD    :    TOA   TAG    32 
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Figure   D.42.       PSD    :    TOA  TAG    48 
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Figure    D.43.       PSD    :    TOA   TAG    51 
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Figure    D.44.       PSD    :    TOA   TAG    64 
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APPENDIX  E 
ADDITIONAL  PLOTS 

These  additional  plots  are  re-f erred  to  in  the  text. 
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Figure  E.l.   GROUP  SPECTRUM  :  TOA  TAG  5 
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Figure  E.2.   GROUP  SPECTRUM  :  TOA  TAG  12 
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Figure  E.3.   GROUP  SPECTRUM  :  TOA  TAG  14 
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Figure  E.4.   GROUP  SPECTRUM  :  TOA  TAG  1? 


224 


>- 


I — 
CO 

a 

UJ 

> 

4L. 

Ul 

a 

on 

UJ 

a. 
en 

on 

UJ 

02 

a 

\n 

a. 

!  en 

csa 

en 

jCM 

1  CM 


I  CM 


CM 
CM 


d^ 


CM 

00 

t— » 

> 

CJ 

UJ 

2 

»— * 

UJ 

"Tf 

ZJ 

*— 4 

a 

UJ 

CM 

CO 

J  Si 

i  *+ 

t 
!co 

CD 
^f 

CM 


th      a 


CD 


co      in 


en      cm 


S3 


CSE90"-!)    3VW  "13d 


Figure  E.5.   PSD  :  TOA  TAG  5  AND  1?  INTERLEAVED. 
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